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ABSTRACT 

•v  ■' 

y 

Inclusion  of  thermal  energy  storage  in  a  pulsed  space 
power  supply  will  reduce  the  mass  of  the  heat  rejection 
system.  In  this  mode,  waste  heat  generated  during  the 
brief  high-power  burst  operation  is  placed  in  the  thermal 
store;  later,  the  heat  in  the  store  is  dissipated  to  space 
via  the  radiator  over  the  much  longer  nonoperational  period 
of  the  orbit.  Thus,  the  radiator  required  is  of  signifi¬ 
cantly  smaller  capacity.  Scoping  analysis  indicates  that 
use  of  lithium  hydride  as  the  thermal  storage  medium  re¬ 
sults  in  system  mass  reduction  benefits  for  burst  periods 
as  long  as  800  s. 
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A  candidate  design  for  the  thermal  energy  storage 
component  utilizes  lithium  hydride  encapsulated  in  either 
304L  stainless  steel  or  molybdenum  in  a  packed-bed  con¬ 
figuration  with  a  lithium  or  sodium-potassium  (NaK)  heat 
transport  fluid. —  Key  issues  associated  with  the  system 
design  include  ptrase-change  induced  stresses  in  the  shell, 
lithium  hydride  aid  shell  compatibility,  lithium  hydride 
dissociation  and  (hydrogen  loss  from  the  system,  void 
presence  and  movement  associated  with  the  melt-freeze 
process,  and  heat  transfer  limitations  on  obtaining  the 
desired  energy  storage  density. 

Preliminary  thermal  analysis  indicates  that,  for 
lithium  hydride  in  a  sphere  of  4-cm  maximum  diameter, 
complete  lithium  hydride  melting  will  be  achieved  in  the 
500  s;  this  provides  a  subunit  with  maximum  energy  storage 
density.  However,  an  elastic  stress  analysis  indicates 
that  the  stresses  generated  in  a  rigid  spherical  shell  by 
the  melting  lithium  hydride  will  necessitate  a  thick  shell, 
significantly  reducing  the  energy  storage  density.  This 
may  be  mitigated  by  cracks  which  form  in  the  lithium 
hydride  during  cooldown,  providing  channels  for  flow  of  the 
melting  salt  during  the  charge  cycle.  Alternative 
approaches  to  shell  stress  reduction  would  control  heat 
flow  so  as  to  locate  the  void  that  forms  in  a  region 
accessible  to  the  melting  lithium  hydride  or  to  use  a 
flexible  shell,  such  as  a  bellows,  for  containment. 

Initial  scoping  experiments  with  encapsulated  lithium 
hydride  have  been  completed  for  cylindrical  cans  of  2.54-cm 
length,  3.81-cm  diameter,  and  0.0635  and  0.0889-cm  (25  and 


2 


35-mil)  wall  thicknesses.  Four  thermal  cycles  were  com¬ 
pleted  without  shell  failure  with  a  0.0889-cm  (35-mil) 
thick  can.  Post-test  examination  of  all  cans  tested  showed 
the  presence  of  numerous  cracks  in  the  lithium  hydride. 


1.  INTRODUCTION 


Power  levels  associated  with  space  defense  weapon  systems  are 
estimated  to  be  very  high,  with  projections  in  the  nultimegawat t  power 
range.  Moreover,  while  these  high  power  levels  are  required  only  in  a 
pulsed  mode,  i.e.,  for  moderately  short  periods  of  time,  the  heat  re¬ 
jection  system  mist  still  be  sized  for  the  full  thermal  load.  Energy 
storage  devices  may  be  employed  to  reduce  the  size  and  mass  of  power 
system  components.  This  report  documents  a  study  on  the  use  of  thermal 
energy  storage  (TES)  in  the  heat  rejection  system  for  this  purpose. 

The  sink-side  TES  unit  accepts  reject  heat  from  the  power 
conversion  system  during  burst  power  operation.  Then,  during  the 
remaining  portion  of  the  orbit,  which  is  typically  an  order  of  magnitude 
longer  than  the  burst  mode  period,  the  stored  reject  heat  is  dissipated 
to  space.  Thus  the  heat  rejection  rate  is  reduced,  and  a  smaller 
radiator  is  required.  In  concepts  where  launch  packages  are  constrained 
by  volume  rather  than  mass,  this  can  be  a  critically  important  attri¬ 
bute.  If  the  reduction  in  radiator  mass  is  greater  than  the  additional 
mass  of  the  storage  module,  overall  system  mass  savings  are  possible. 

Lithium  hydride,  because  of  its  superior  heat  storage  properties 
and  convenient  melting  temperature,  appears  to  be  the  best  candidate  for 
use  in  sink-side  TES  applications.  To  maximize  storage  density,  both 
sensible  and  latent  heat  characteristics  of  the  medium  are  used. 
Several  TES  design  concepts  based  on  LiH  storage  are  being  evaluated 
analytically  and/or  experimentally.  One  concept  employs  the  direct 
contact  of  the  heat  transport  medium  (NaK  or  lithium)  with  lithium 
hydride  to  effect  the  highest  heat  transfer  rates  and  storage  densities 
possible.1  Another  approach  uses  encapsulated  shapes  of  lithium  hydride 
in  a  packed  bed  or  bulk  lithium  hydride  can  be  used  in  a  conventional 
shell  and  tube  configuration. 
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This  report  focuses  on  the  encapsulated  shapes  concept  and 
describes  the  analytical  and  experimental  work  that  has  been  completed 
to  date.  Since  this  program  is  in  its  early  stages,  most  of  the  work 
has  been  scoping  in  nature.  However,  some  fairly  detailed  phase-change 
heat  transfer  analysis  has  been  performed  in  spherical  geometry.  In 
addition,  analytical  efforts  have  been  extended  to  characterizing  the 
behavior  of  the  void  which  occurs  due  to  the  large  lithium  hydride 
volumetric  shrinkage  during  the  solidification  process. 

Analytical  and  experimental  work  have  addressed  basic  feasibility 
issues  associated  with  encapsulated  lithium  hydride:  phase-change 
induced  stresses  on  heat-up,  hydrogen  diffusion  through  shell  walls, 
heat  transfer  rate,  void  behavior,  and  material  constraints.  The  impact 
of  these  issues  on  the  design  of  lithium  hydride  encapsulated  shapes  has 
been  evaluated,  and  design  alternatives  have  been  identified  for  circum¬ 
venting  the  key  problem  areas. 

The  following  sections  of  this  report  contain  a  description  of  the 
encapsulated  shape  design  concept  a  presentation  of  the  analytical  wcrk 
completed,  results  of  an  assessment  of  the  value  (i.e.,  power  system 
mass  and  area  savings)  of  sink-side  thermal  energy  storage  are  pre¬ 
sented,  and  an  evaluation  of  the  feasibility  issues.  This  is  followed 
by  a  discussion  of  the  initial  scoping  experiments  with  encapsulated 
lithium  hydride.  Finally,  based  on  work  to  date,  preliminary  con¬ 
clusions  are  provided. 


2.  SYSTEM  CONCEPT  ANALYSIS 


A  preliminary  system  value  analysis  was  performed  to  determine  if 
the  heat  storage/rejection  concept  had  sufficient  merit  to  warrant 
development.  The  analysis  also  was  useful  in  determining,  in  a  pre¬ 
liminary  manner,  the  minimum  required  storage  density  that  resulted  in 
the  encapsulated  lithium  hydride  storage  unit  concept  being  feasible  for 
burst  power  applications.  A  scoping  design  analysis  was  then  conducted 
to  determine  if  the  storage  unit  could  be  designed  to  meet  the  storage 


density  goals. 


2. 1  System  Concept 


A  simplified  schematic  of  a  no  effluent  space  power  system 
employing  encapsulated  sink-side  thermal  storage  is  shown  in  Fig.  2.1. 
Depending  upon  the  application  the  heat  source  could  be  nuclear,  chem¬ 
ical  or  solar.  Power  conversion  can  be  accomplished  using  dynamic 
systems  such  as  Rankine  and  Brayton  cycles,  or  static  systems  such  as 
thermionics.  As  mentioned  previously  waste  heat  would  be  stored  in  the 
encapsulated  lithium  hydride  during  burst  power  operation.  The  stored 
heat  would  then  be  rejected  over  the  longer,  non-operational  portion  of 
the  orbit.  The  candidate  heat  transport  fluids  include  NaK  and  lithium, 
while  the  radiator  could  utilize  heat  pipe  technology.  Other  tech¬ 
nologies  as  they  are  developed  would  provide  alternatives. 

ORNL-  DWG87  4004  ETO 
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Fig.  2.1.  Schematic  of  power  system  employing  encapsulated  sink- 
side  thermal  storage. 


The  mass  of  the  heat  rejection  system  using  storage  was  compared  to 
the  mass  of  the  baseline  system,  which  used  only  a  radiator,  to  give  a 
measure  of  the  value  of  the  storage  system.  In  the  baseline  system, 
waste  heat  is  rejected  immediately  via  radiation  during  the  power  gener¬ 
ation  portion  of  the  orbit.  Thus,  the  radiator  is  sized  to  handle  the 
instantaneous  load  imposed  by  the  power  cycle.  With  thermal  storage  the 
reject  heat  is  placed  into  storage  during  the  power  generation  portion 
of  the  orbit;  then,  during  the  non-operational  portion  of  the  orbit  the 
stored  heat  is  rejected  via  radiation. 

The  major  parameters  of  concern  include:  energy  storage  density, 
minimum  temperature  of  the  storage  medium,  burst  operation  period  and 
radiator  specific  mass.  The  storage  density  depends  upon  the  design  of 
the  unit  (which  is  addressed  in  the  subsequent  section)  and  the  tempera¬ 
ture  swing  (the  difference  between  the  maximum  and  minimum  storage  tem¬ 
perature).  The  storage  medium  utilizes  the  sensible  heat  of  the  liquid 
and  solid  phases  in  addition  to  the  latent  heat  of  fusion.  Thus,  the 
total  heat  stored  is  strongly  influenced  by  the  temperature  swing.  The 
maximum  temperature  of  the  storage  unit  is  controlled  by  the  reactor 
coolant  temperature,  which  in  this  analysis  was  assumed  to  be  1100  K. 
The  temperature  swing,  therefore,  was  fixed  by  the  minimum  storage  tem¬ 
perature.  Variation  of  the  minimum  storage  temperature  produces  two 
opposing  effects.  As  the  minimum  storage  temperature  increases,  the 
temperature  swing  decreases.  This  decreases  the  storage  density  and  in¬ 
creases  the  required  storage  mass.  However,  the  average  radiator  flux 
increases  and  this  decreases  the  required  radiator  area  and  mass. 

The  total  orbit  time  was  fixed  at  6000  s  and  the  burst  period 
varied  between  100  and  900  s.  Thus,  the  non-operational  period  (the 
time  available  for  heat  rejection)  varied  between  5900  and  5100  s.  The 
baseline  (i.e.,  non-storage)  system  radiator  was  assumed  to  operate  at 
1100  K.  This  results  in  a  radiator  heat  flux  of  74.)  kW/m2.  The 
radiator  specific  mass  was  fixed  at  20  kg/m2,  which  is  felt  to  be 
representative  of  a  survlvable  heat  pipe  system. 


Previous  analysis  by  Olszewski  and  Morris4  indicated  that  the 
optimal  heat  storage/rejection  system  operated  with  minimum  storage 
temperature  in  the  range  of  500  to  700  K.  The  analysis  was,  therefore, 
confined  to  this  range  of  lower  operating  temperatures. 

The  value  analysis  results  are  given  in  Fig.  2.2.  The  crossover 
time  is  defined  as  the  time  at  which  the  storage  heat  rejection  and  the 
radiator-only  systems  are  of  equal  mass.  For  total  generation  times 
less  than  the  crossover  time  the  storage  system  is  lighter.  Thus,  a 
storage  heat  rejection  system  with  a  storage  density  of  3  MJ/kg  operat¬ 
ing  with  a  lower  temperature  of  700  K  will  show  advantage  over  a 
radiator-only  system  for  generation  times  of  600  s  or  less.  \s  antici¬ 
pated  but  not  illustrated  in  the  figure  (see  Ref.  2  for  details)  the 
value  of  storage  (i.e.,  mass  savings)  increases  with  decreasing  genera¬ 
tion  time.  To  have  reasonable  applicability  to  burst  power  needs  it  is 
desirable  that  the  storage  system  have  a  crossover  time  of  at  least 
500  s.  From  Fig.  2.2  it  is  evident  that  this  will  require  a  system 
storage  density  of  at  least  3  MJ/kg.  Thus,  the  research  program  is 
directed  at  achieving  a  system  storage  density  exceeding  this  value. 
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Fig.  2.2.  Heat  rejection  system  crossover  time  for  varying  minimum 
storage  temperatures. 


2.3  Storage  System  Conceptual  design 


Preliminary  system  conceptual  designs  were  prepared  for  the  pro¬ 
posed  storage  system  using  encapsulated  lithium  hydride  shapes  in  a 


packed  bed  configuration.  The  purpose  of  this  design  analysis  was  to 
determine  if  the  system  storage  densities  could  meet  the  goal  of  3  MJ/kg 
or  greater  as  determined  by  the  value  analysis.  Conceptual  designs  were 
completed  for  modules  that  would  meet  the  requirements  of  the  following 
mission: 

•  burst  generation  time,  500  s 

•  system  output  power,  300  MW(e) 

•  power  system  efficiency,  25% 

•  reject  heat  load,  900  MW(t) 

•  heat  stored,  450  GJ. 

The  TES  system  operational  storage  density  depends  upon  many  design 
related  factors.  The  storage  medium  used,  thickness  of  the  encapsulat¬ 
ing  shell,  packing  density  of  the  encapsulated  shapes  and  choice  of 
coolant  have  significant  impact  on  the  mass  requirements  of  the  non¬ 
storage  components  of  the  unit  and  this  affects  the  storage  density  of 
the  system. 

System  storage  densities  were  calculated  for  several  design  op¬ 
tions.  Lithium  and  NaK  were  examined  as  the  heat  transport  fluids. 
Also,  the  storage  medium  was  examined  using  naturally  occurring  lithium 
(Li7)  and  the  isotope  lithium  6  (Li6).  Preliminary  information  in¬ 
dicates  that  Li6  possess  the  same  molar  thermal  properties  as  Li7. 
Thus,  on  a  mass  basis  the  heat  of  fusion  and  specific  heat  will  be  12% 
higher  for  the  Li6. 

It  was  assumed  that  the  lithium  hydride  was  formed  in  spheres  and 
encapsulated  with  a  0.013-cm  (5-mil)  thick  stainless  steel  shell  (this 
shell  thickness  assumption  represents  a  design  goal).  As  described  in 
Sect.  3.2,  lithium  hydride  heat  transport  properties  limit  the  maximum 
lithium  hydride  sphere  diameter  to  3.8  cm  (the  maximum  sphere  diameter 
that  allows  for  complete  melting  of  the  lithium  hydride  in  the  500  s 
burst  time).  With  this  design  the  mass  of  the  lithium  hydride  sphere 
(using  Li7)  is  15.8  g  and  that  of  the  steel  shell  is  4.7  g. 

Results  of  the  design  analysis  are  presented  in  Table  2.1.  As  in¬ 
dicated  by  these  results  lithium  is  the  preferred  thermal  transport 
medium  since  all  designs  yield  storage  system  densities  in  excess  of 


5:1 


$ 


03 

2  ^ 

d  a> 

S  o-  ^ 

t-»  e  ^ 

c  <u  w 


r- 

*— H 

ro 

00 

co 

o 

vO 

CN 

r-. 

* 

• 

• 

• 

• 

LTl 

m 

<j- 

<r 

ro 

p*- 

CN 

CO 

O' 

vD 

<r 

oc 

• 

• 

• 

• 

• 

co 

rn 

co 

rn 

CN 

vD 

O 

m 

CN 

vO 

r>. 

<r 

O 

sC 

• 

• 

• 

• 

• 

in 

m 

in 

<r 

<r 

CO 

O 

, 

CO 

, 

r^ 

m 

CN 

00 

m 

• 

• 

• 

• 

• 

<r 

<r 

<r 

rn 

CO 

CN 

00 

<T 

vT> 

CO 

o 

vC 

m 

O' 

in 

• 

• 

• 

• 

• 

in 

<r 

<T 

C' 

CO 

o 

n» 

co 

r-v 

00 

r^ 

CN 

O*' 

• 

• 

• 

« 

• 

co 

m 

co 

CN 

CN 

CM 

o> 

in 

r- 

<r 

co 

CTN 

vD 

CN 

00 

• 

• 

• 

• 

• 

m 

<3- 

<■ 

CO 

i o 

O 

CN 

CN 

<T 

CN 

O' 

\-D 

m 

• 

• 

• 

• 

• 

<t 

<5* 

CO 

CO 

CO 

<T 

*”H 

v£> 

<r 

CO 

O 

m 

o> 

co 

v£> 

• 

• 

• 

• 

• 

O' 

oo 

n- 

r-- 

MO 

o 

<r 

vO 

CN 

O 

cr 

<J- 

O' 

<r 

oo 

r^. 

N 

lO 

vO 

m 

• 

• 

• 

o 

o 

o 

o 

O 

o 

o 

o 

o 

O 

co 

<r 

m 

o 

r~» 

-♦*.-**>* 


asK 


3  MJ/kg.  The  use  of  Li6  is  also  preferred  because  of  Its  enhanced  stor¬ 
age  density.  As  indicated  by  the  detailed  mass  breakdown  information 
presented  in  Table  2.2  the  mass  of  the  encapsulated  storage  medium 
(lithium  hydride  and  shell)  and  the  coolant  essentially  determine  the 
system  mass.  Since  the  storage  capacity  of  the  coolant  is  much  less 
than  that  of  the  encapsulated  lithium  hydride  (<10%  of  the  storage 
capacity  of  the  lithium  hydride),  increasing  the  packing  density  (by 
using  two  sphere  sizes)  has  a  significant  impact. 


Table  2.2.  Storage  system  detailed  mass 
breakdown  for  700  K  minimum  storage 
temperature  design  option 


Storage  medium 
lithium  isotope 

Coolant 

Packing 

densitv 

(%) 

Mass 

(kg) 

St  orage 
medium 

Coolant 

Pi  ping 
con t  a i ne  r 

etc. 

Total 

Li7 

NaK 

60 

95,177 

68,316 

4,500 

168,194 

Li6 

NaK 

60 

86,225 

68,471 

4,500 

159,196 

Li  ° 

Li 

60 

79,435 

35,146 

4,500 

1  19,081 

Li6 

Li 

75 

84,870 

18,776 

4,500 

108 ,146 

Using  the  results  in  Table  2.1  with  those  in  Fig.  2.2  confirm  the 
applicability  of  the  packed  bed  design  for  burst  power  needs.  The 
design  using  Li^,  a  packing  density  of  75%  and  lithium  coolant  yields  an 
operational  storage  density  of  4.2  MJ/kg  at  a  minimum  storage  tempera¬ 
ture  of  700  K.  Combining  these  values  with  the  value  analysis  results 
in  Fig.  2.2  indicates  that  the  storage  system  would  be  attractive  for 
burst  times  as  long  as  800  s.  This  greatly  exceeds  the  original  goal  of 
500  s. 

In  summary,  system  storage  densities  of  3  MJ/kg  or  greater  will  be 
required  to  produce  a  feasible  sink-side  heat  storage  unit.  It  is  rea¬ 
sonable  to  expect  operational  storage  densities  in  excess  of  3  MJ/kg 
using  the  encapsulated  lithium  hydride  packed  bed  concept.  The  most 
promising  system  uses  Li6  as  the  storage  medium  with  a  packing  density 
of  75%  and  lithium  as  the  thermal  transport  medium. 
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3.  DEVELOPMENT  ISSUES  AND  EVALUATION 


3.1  Overview  of  Issues 

As  indicated  in  Sect.  2.3,  the  lithium  hydride  encapsulating  shell 
must  be  sufficiently  thin  to  achieve  a  storage  system  design  which  pro¬ 
vides  overall  system  mass  savings.  However,  there  are  several  factors 
(design  issues)  which  might  not  allow  the  use  of  a  thin  shell,  and 
therefore  need  to  be  investigated.  The  key  issue  of  concern  is  stress 
induced  in  the  shell  during  heatup.  The  stress  results  from  hydrostatic 
forces  as  lithium  hydride  expands  during  phase  change.  A  second  area  of 
concern  is  hydrogen  diffusion  through  the  shell.  The  rate  of  hydrogen 
loss  varies  inversely  with  shell  thickness.  Thus,  hydrogen  loss 
increases  with  decreasing  shell  thickness. 

As  mentioned  previously,  if  lithium  hydride  is  encapsulated  in 
spheres,  the  poor  heat  transport  properties  of  lithium  hydride  limit  the 
maximum  sphere  size  to  about  3.8  cm  (for  a  500  s  burst  time).  Void 
behavior  plays  a  key  role  in  assessing  the  design  limitations  resulting 
from  poor  lithium  hydride  heat  transport,  and  even  more  importantly,  on 
phase-change  induced  shell  stress.  A  void  at  a  shell  surface  reduces 
heat  transfer  although  hydrogen  is  a  good  gas  thermal  conductor.  The 
way  in  which  the  void  influences  shell  stress  is  described  later  in  this 
section. 

From  a  materials  standpoint,  compatibility  of  the  shell  with 
lithium  hydride,  lithium,  and  hydrogen  is  required  over  a  fairly  wide 
range  of  temperatures.  However,  compatibility  concerns  are  somewhat 
mitigated  by  the  rather  short  projected  system  integrated  operating 
time.  Other  materials  concerns  include  thermal  shock  resistance, 
material  cost,  ease  of  fabrication,  strength,  and  ductility. 

Sections  3.2  through  3.6  describe  models,  analysis,  and  research 
conducted  to  date  to  assess  and  resolve  the  key  design  issues. 
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3.2  Thermal  Analysis 


3.2.1  Description  of  heat  transfer  model 

A  two-dimensional  finite  difference  heat  transfer  computer  code  is 
being  developed  to  model  the  thermal  performance  of  lithium  hydride  en¬ 
capsulated  in  a  spherical  shell.  The  model  uses  the  "enthalpy"  method 
to  account  for  phase  change  in  the  salt.  In  the  enthalpy  method,  the 
location  of  the  phase  change  front  is  not  determined  directly  but  is  in¬ 
ferred  from  the  energy  content  of  the  material.  This  method  has  been 
used  extensively  at  ORNL,  for  example  see  Refs.  3  and  4.  A  non-typical 
aspect  of  the  present  computer  program  is  that  the  fully  implicit  form 
of  the  heat  balance  equations  is  solved,  rather  than  the  explicit  form 
usually  employed.  The  two  main  reasons  for  using  the  fully  implicit 
solution  are  to  avoid  stability  problems  associated  with  the  explicit: 
solution  and  to  minimize  the  difficulty  in  implementing  a  natural  con¬ 
vection  model  during  the  latter  part  of  the  project.  Heat  balance  equa¬ 
tions  are  solved  on  a  regular  grid  in  r,  9  coordinates.  The  equations 
are  solved  directly  for  the  one-dimensional  case  or  by  line  relaxation 
for  the  two-dimensional  case.  The  solution  method  is  given  in  Ref.  5. 
Temperatures  are  scaled  so  the  scaled  temperature,  T' ,  at  the  melt  point 
Is  zero: 

T'  =  T  —  T 

mp 

where 

T  -  956  K  (1260°F) 

mp 

Energy  content  is  scaled  so  the  energy  content,  E,  of  the  solid  at  the 
melt  point  is  zero.  Therefore, 


E  -  C  T 
s 


E  -  V'  +  Hs* 


For  T'  <  0 


For  T'  >  0 


where  Cg  and  are  the  specific  heats  of  the  solid  and  liquid,  and  Hg^ 
is  the  heat  of  fusion  [“  2,6  MJ/kg  (“1109  Btu/lbm)].  When  a  given  node 


As  can  be  seen  from  the  plot,  agreement  between  the  two  solutions 
was  satisfactory.  Results  of  two-dimensional  calculations  were  checked 
against  results  from  the  explicit  HEATING6  code  (see  Ref.  A)  and  were 
found  to  agree  within  1%,  thus  confirming  that  the  numerical  formulation 
was  correct.  Comparison  with  experimental  results  will  be  needed  to 
validate  the  model. 


3.2.2  Heat  transfer  analysis 

The  finite  difference  code  was  used  to  determine  the  thermal  re¬ 
sponse  of  single  spheres  of  various  diameters  exposed  to  liquid  metal 
convection  when  the  fluid  temperature  undergoes  a  step  change  from  300  K 
(80°F)  to  1100  K  (1520°F).  This  simulates  the  power-up  of  an  orbital 
power  system  where  spheres  containing  lithium  hydride  arranged  in  a 
packed  bed  would  be  used  as  a  heat  sink  with  NaK  or  lithium  as  the  heat 
transfer  fluid.  Theoretical  relationships  for  convection  coefficients 
for  a  liquid  metal  flowing  past  a  single  sphere  have  been  developed. 
The  convection  coefficient  for  a  single  sphere  varies  strongly  with 
angle  from  the  leading  stagnation  point,  with  a  maximum  at  the  leading 
stagnation  point  and  a  minimum  at  the  trailing  stagnation  point.  How¬ 
ever,  it  was  felt  that  an  average  value  over  the  entire  surface  would  be 
more  representative  of  convection  in  a  packed  bed.  Based  on  correla¬ 
tions  given  in  Ref.  7,  an  average  value  of  17  kW/m2-K  (3000  Btu/h-ft2- 
°F)  [for  a  fluid  velocity  of  0.30  m/s  (1  ft/s)]  was  used  in  the  study. 

The  mass  of  lithium  hydride  in  the  spheres  was  set  so  they  are  com¬ 
pletely  filled  when  the  lithium  hydride  is  liquid.  The  shells  are 
stainless  steel,  with  thickness  equal  to  1%  of  the  inside  radius  [e.g. , 
for  R^  *  1.9  cm,  shell  thickness  -  0.019  cm  (7.5  mil)].  Linear  scaling 
of  shell  thickness  can  be  justified  due  to  stress  considerations,  and  as 
a  result  the  shell  has  no  effect  when  comparing  the  specific  energy  con¬ 
tent  of  different  size  spheres.  Properties  for  the  shell  material,  and 
solid  and  liquid  lithium  hydride  were  considered  to  be  constant.  A  list 
of  property  values  used  in  the  study  is  given  in  Table  3.1. 

In  order  to  formulate  a  thermal  model  for  melting  of  lithium 
hydride  in  a  sphere,  two  important  assumptions  must  be  made.  The  first 
of  these  is  the  initial  location  of  the  void  in  the  solid  material.  The 
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Table  3.1.  Property  values  used  in  thermal  analysis 


Material 

Density,  g/cm3 
(lbm/ft3) 

Specific  heat 
J/g-K 

(Btu/lbm-°F) 

Thermal  conductivity 
J/s -cm-K 
(Btu/h-f t-°F) 

Solid  LiH 

0.69  (43.08) 

6.28  (1.50) 

0.042  (2.42) 

Liquid  LiH 

0.55  (34.30) 

7.37  (1.76) 

0.021  (1.21) 

Stainless  steel 

7.81  (487.6) 

0.582  (0.139) 

0.22  (12.86) 

LiH  melt  temperature  =  956  K  (1260°F). 

LiH  heat  of  fusion  =  2.58  MJ/kg  (1109  Btu/lbm). 

behavior  of  the  void  during  solidification  is  not  well  understood,  and 
developing  an  understanding  of  void  mechanics  is  one  of  the  main  goals 
of  this  project  (void  mechanics  is  discussed  in  detail  in  Sect.  3.3). 
Minimum  energy  considerations  indicate  that  in  the  case  of  uniform 
cooling  in  micro-gravity,  a  spherical  void  will  form  in  the  center  of 
the  container  (no  stability  analysis  has  been  done  to  determine  if  other 
shapes  are  possible).  This  assumption  is  made  in  the  present  study.  A 
second  critical  assumption  that  must  be  made  concerns  how  the  volumetric 
expansion  due  to  melting  is  accounted  for  in  the  case  where  the  liquid 
is  not  in  contact  with  the  void.  The  first  possibility  is  that  the 
shell  expands  to  provide  the  required  volume.  This  case  is  not  of 
Interest  in  the  heat  transfer  study,  since  any  practical  shell  would  be 
able  to  expand  very  little  before  rupturing.  It  is,  however,  of  great 
interest  in  determining  the  design  of  the  shell.  This  problem  is 
discussed  in  detail  in  Sect.  3.4.4.  A  second  way  the  volume  expansion 
could  be  accommodated  is  for  the  solid  to  be  crushed  uniformly  inward 
into  the  void.  Whether  or  not  this  is  a  realistic  possibility  with  a 
practical  shell  thickness  depends  on  the  compressive  strength  of  the 
solid  lithium  hydride,  but  in  any  event  it  provides  a  limiting  case  for 
the  heat  transfer  study.  A  third  possibility  is  for  the  excess  liquid 
to  leak  through  the  solid  into  the  void.  This  would  likely  occur  if 
there  are  sufficient  cracks  in  the  solid.  Heat  transfer  calculations 
were  carried  out  assuming  the  second  case  mentioned  above,  hereafter 
referred  to  as  the  "crush"  model,  and  the  last  case,  hereafter  referred 
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to  as  the  "leak"  model.  In  the  leak  model  calculations,  the  means  by 
which  the  liquid  makes  its  way  through  the  solid  is  not  considered.  It 
is  assumed  that  liquid  at  the  fusion  temperature  appears  in  a  freezing 
volume  at  the  inner  surface  of  the  solid.  In  reality,  the  phase  front 
at  the  inner  surface  of  the  solid  may  not  be  uniform  since  the  liquid 
passes  through  localized  cracks  in  the  solid.  The  processes  involved  in 
the  two  models  are  illustrated  in  Fig.  3.2. 

Figure  3.3  compares  the  fraction  of  the  lithium  hydride  melted 
versus  time  for  the  crush  and  the  leak  models.  The  results  shown  are 
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Fig.  3.2.  Crush  and  leak  thermal  models. 

a.  Crush  model  —solid  is  uniformly  crushed  toward  center  to  account 
for  volume  expansion  at  melt  front. 

b.  Leak  model  —excess  liquid  from  melting  volume  appears  in  freezing 
volume  at  inner  surface  of  solid.  How  it  gets  there  is  not 
considered. 
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Fig.  3.3.  Fraction  melted  versus  time,  RQ  “  1.9  cm. 


for  a  1.9-cm  radius  sphere.  During  the  early  stage  of  the  melting 
process,  the  two  models  give  very  similar  results,  but  as  the  process 
continues  the  rate  of  melting  predicted  by  the  leak  model  begins  to 
exceed  that  of  the  crush  model.  This  difference  increases  with  time. 
This  result  is  due  primarily  to  two  effects.  The  conduction  path 
through  the  liquid  is  shorter  in  the  leak  model.  In  addition,  some  heat 
is  carried  along  with  the  liquid  that  leaks  through  the  solid  and  thus 
bypasses  the  conduction  path  through  the  solid. 

Figure  3.4  shows  a  comparison  of  energy  gain  vs  time  for  spheres 
ranging  from  0.635  to  3.81  cm  in  radius.  Assuming  a  total  time  of  500  s 
available  for  charging,  the  1.9-cm  radius  sphere  Is  the  largest  that 
attains  its  maximum  potential  energy  storage  In  the  required  time.  As 
shown  in  Fig.  3.5,  the  3.81-cm  radius  sphere  requires  about  800  s  to 
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Fig.  3.4.  Energy  storage  versus  time  for  lithium  hydride  spheres 
(leak  model). 


completely  melt.  Table  3.2  lists  energy  content  after  500  s  for  the 
various  sphere  sizes  predicted  by  the  crush  and  leak  models.  It  is 
interesting  to  note  that  the  results  for  the  two  cases  do  not  differ 
greatly  (the  maximum  difference  is  about  2%). 

In  summary,  the  mechanism  by  which  the  volume  expansion  due  to 
melting  is  accommodated  is  not  critical  in  the  heat  transfer  model,  as 


Table  3.2.  Energy  content  at  500  s 
for  various  sphere  sizes 


Energy  content  at  t  =  500  s 


Diameter 


Crush  model 


Leak,  model 


(inches)  (MJ/kg)  ("  max)  (MJ/kg)  (%  max) 


0.0  5.559 
9.96  5.558 
9.50  5.550 
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Fig.  3.5.  Fraction  melted  versus  time  for  lithium  hydride  spheres 
(leak  model). 

the  crush  and  leak  models  give  very  similar  results.  Both  models 
indicate  that  for  complete  melting  in  500  s  with  an  1100  K  source 
temperature,  the  maximum  container  diameter  that  can  be  used  is  about 
3.81  cm. 


3.3  Void  Mechanics 

3.3.1  Physics  of  solid-llquld-vapor  systems 

Determination  of  liquid-vapor  interface  shapes  is  a  class  of 
problems  which  has  received  a  great  deal  of  study.  However,  with  the 
exception  of  some  relatively  simple  geometries,  no  straightforward 
method  of  solving  for  interface  shapes  is  available.  Surface  energy 
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arises  because  molecules  at  the  surface  of  a  liquid  are  not  attracted 
equally  from  all  sides  as  are  those  in  the  bulk  liquid.  Thus,  these 
molecules  have  some  excess  energy  compared  to  molecules  in  the  bulk 
liquid.  The  unequal  attraction  experienced  by  the  surface  molecules  re¬ 
sults  in  a  force  acting  along  the  liquid  surface.  This  force  can  be 
represented  as  a  uniform  tension  acting  on  a  hypothetical  membrane 
stretched  over  the  surface  and  is  referred  to  as  surface  tension,  a. 
The  relationship  of  the  surface  tensions  between  the  liquid  and 
vapor,  a£v»  the  solid  and  vapor,  a3w>  and  the  solid  and  liquid,  o  is 
expressed  by  the  equation 


o  —  a  „ 
sv  si 
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justified  on  that  basis,  however,  since  the  force  normal  to  the  solid 
surface  is  not  balanced.  Alternately,  Young's  equation  can  be  derived 
based  on  thermodynamic  considerations.8 

Although  the  surface  tension  and  contact  angle  are  manifestations 
of  microscopic  forces,  together  they  have  been  found  to  provide  an 
adequate  representation  of  the  macroscopic  behavior  of  capillary  sys¬ 
tems.  The  main  disadvantage  of  the  surface  tension-contact  angle 
representation  is  that  the  contact  angle  can  be  strongly  influenced  by 
the  condition  of  the  surface  and  is  subject  to  hysteresis  when  the 
liquid  advances  or  recedes,  thus  causing  some  uncertainty  in  calcula- 
t ions . 

There  are  two  basic  approaches  that  can  be  used  to  calculate  inter¬ 
face  shapes.  The  first  involves  minimizing  the  free  energy,  F,  of  the 
system,  which  for  an  isothermal  capillary  system  of  mass,  m,  is  the  sum 
of  the  surface  energy  and  potential  energy,  and  can  be  expressed  as: 

F  =  JA  a  dA  +  /v  pgh  dV 
c 

where  Ac  is  the  capillary  area,  V  is  the  total  system  volume  and  h  is 
the  elevation  relative  to  a  fixed  reference.  The  equation  for  the  mini¬ 
mum  energy  surface  can  be  determined  by  applying  the  calculus  of  varia¬ 
tions  (<5F  =  0)  with  the  constraint  that  the  liquid  volume  is  constant. 
This  method  has  been  applied  successfully  to  determine  equilibrium  in¬ 
terface  shapes  in  axisymmetric  containers.8 

A  second  method  that  can  be  used  to  calculate  interface  shapes  in¬ 
volves  balancing  forces  on  the  interface.  Consider  a  small  surface  ele¬ 
ment  dA  of  the  liquid-vapor  interface  as  shown  in  Fig.  3.7.  Element  dA 
has  minimum  and  maximum  radii  of  curvature  Rj  and  R2  with  origins  along 
the  normal  to  dA.  Since  they  represent  principal  (maximum  and  minimum) 
curvatures,  the  planes  of  R^  and  R2  are  perpendicular  to  each  other. 
The  length,  L^,  of  dA  in  the  plane  of  Rj  is  R^A0^,  where  '  s  assumed 
to  be  very  small.  Similarly,  the  length  L2  of  dA  in  the  plane  of  R2 
is  R2&02.  In  order  for  forces  directed  along  the  normal  to  be  balanced, 
the  normal  component  of  the  surface  tension  forces  acting  on  the  edges 
of  dA  must  equal  the  pressure  force  on  dA,  The  equation  expressing  the 


sin  A0  »  A0 ,  the  equation  becomes 


A02 

(Pl  -  P2)(R1A61)(R2A02)  =  ZCRjAQ^  a  —  +  2  (R2A62) 


A0j 

a - 

2 


which  reduces  to 


PX  -  P2  =  o 


(3.1) 


Equation  (3.1)  is  called  the  Laplace  equation  and  is  the  classical 
expression  for  equilibrium  of  a  liquid-vapor  interface.  It  should  be 
noted  that  the  Laplace  equation  can  also  be  derived  from  minimum  energy 
considerations.10  Of  course,  both  of  the  methods  mentioned  above  must 
give  the  same  result,  so  the  choice  of  method  depends  on  the  particular 
problem  being  considered.  Most  sources  indicate  that  the  minimum  energy 
formulation  is  more  straightforward  for  most  situations,  with  the  excep¬ 
tion  of  a  few  simple  cases.  However,  it  appears  that  the  number  of 
cases,  simple  or  otherwise,  for  which  a  closed  form  solution  can  be 
obtained  Is  very  limited. 


3.3.2  Determination  of  interface  shape 
for  an  axisymmetric  case 

Fortunately,  for  an  axisymmetric  case  a  vertical  force  balance  can 
be  used  to  find  the  interface  shape.  Consider  an  element  of  the  inter¬ 
face  with  length  dL,  as  shown  in  Fig.  3.8.  This  element  forms  an 
annular  ring  about  the  axis  of  symmetry.  The  pressure  difference  across 
the  interface  is  Pv  —  P^,  where  Py  and  P^  are  the  vapor  and  liquid 
pressures.  P^  can  be  expressed  in  terms  of  the  liquid  pressure  at  the 
origin,  P^q ,  and  the  hydrostatic  head;  so 

P  —  P  =  P  —  P  +  pgz 
vZvlo 

The  projected  area  that  the  pressure  difference  acts  over  is  27TrAr. 
Therefore  the  pressure  force  in  the  -z  direction  is 


(P  -  P  +  pgz)  2ti  rAr 


24 


25 


Eqs.  (3.5)  and  (3.6)  can  be  solved  by  the  Runge-Kutta  method  to  yield 
Interface  shapes  for  given  values  of  o  ,  P  and  g.  The  solution  is 
iterative,  with  values  of  RQ  being  guessed  to  provide  the  required  vapor 
volume  and  contact  angle  between  the  liquid  and  container.  Equations 
similar  to  Eqs.  (3.5)  and  (3.6)  were  used  in  Ref.  11  to  develop  an 
extensive  set  of  curves  that  give  meniscus  shape  as  a  function  of 
Pg/°£v  an<i  Figure  3.9  shows  interface  shapes  for  liquid  lithium 

hydride  in  a  spherical  container  for  void  fractions  ranging  from  5  to 
25%.  It  should  be  noted  that  for  the  micro-gravity  case  the  interface 
shape  is  a  segment  of  a  sphere,  and  for  the  case  where  0=0°,  shown  in 
Fig.  3.9c,  the  Interface  is  a  spherical  bubble  which  has  no  affinity  for 
the  container  wall  and  may  be  completely  surrounded  by  the  liquid. 

3.3.3  Effect  of  surface  tension  gradients 

In  order  to  determine  the  shape  and  location  of  the  void  in  the 
above  case,  it  is  necessary  to  consider  the  effect  of  surface  tension 
gradients.  Surface  tension  gradients  can  be  caused  by  variations  in  the 
chemical  composition  of  the  liquid,  and  surface  tension  in  general 
decreases  linearly  with  increasing  temperature.  A  surface  tension 
gradient  produces  a  shear  stress  along  the  surface  of  the  liquid,  caus¬ 
ing  liquid  to  move  toward  the  region  of  higher  surface  tension.  Conse¬ 
quently,  the  bubble  tends  to  move  toward  the  region  of  lower  surface 
tension.  This  phenomenon  can  be  explained  in  terms  of  the  surface 
energy  as  well,  since  the  surface  energy  is  decreased  when  the  surface 
temperature  is  higher.  Flow  resulting  from  surface  tension  gradients  is 
called  Marangoni  flow.  The  complete  vector  equation  describing 
Marangoni  flow,  given  in  Ref.  10,  is 


where  n  is  the  unit  normal  to  the  surface  and  T  is  the  viscous  stress 
tensor  at  the  surface  of  the  liquid,  not  including  pressure  terms. 
Although  the  above  equation  would  be  very  difficult  to  apply  for  a 
general  case,  it  shows  that  if  a  surface  tension  gradient  exists  there 
must  be  motion  in  the  fluid. 
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The  effect  of  a  surface  tension  gradient  on  the  interface  shape  can 
be  estimated  using  Eq.  (3.2).  In  the  present  case,  the  surface  tension 
will  vary  linearly  in  the  z  direction,  so 


O.  -  0  +  \z 

Iv  o 


where  Oq  is  the  surface  tension  at  the  origin  and  X 


■  — .  Inserting 


this  into  Eq.  (3.2)  yields 
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Following  the  same  stepB  used  to  develop  Eq.  (3.6)  the  resulting 
differential  equation  is 
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(3.7) 


Eq.  (3.7)  can  be  substituted  for  Eq.  (3.6)  and  used  together  with 
Eq.  (3.5)  to  solve  for  the  interface  shape.  Figure  3.10  shows  a  com¬ 
parison  of  the  predicted  void  shape  for  a  zero  g,  zero  contact  angle 
situation  with  no  surface  tension  gradient  versus  the  predicted  shape 
with  a  surface  tension  gradient  corresponding  to  a  22  K/cm  (100°F/in. ) 
temperature  increase  in  the  -z  direction.  As  shown  in  the  figure,  the 
effect  of  the  surface  tension  gradient  on  the  interface  shape  is  very 
small,  resulting  in  a  slight  stretching  of  the  bubble  along  the  z-axis. 

The  effect  of  the  surface  tension  gradient  on  the  bubble  location 

is  much  more  significant.  An  estimate  of  the  resultant  force  on  a 

spherical  bubble  due  to  Marangoni  flow  is  developed  in  Ref.  12  for  the 

case  of  a  uniform  surface  tension  gradient.  The  shear  stress  due  to  the 

surface  tension  gradient  is  integrated  over  the  surface  of  the  bubble, 

and  the  resultant  force,  F  .  in  the  z  direction  is  found  to  be 

nr 
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where  RQ  Is  the  bubble  radius.  The  buoyancy  force,  F^,  is  given  by 


Fb  *  1  "V  (p£  ~  \)  g 

Under  normal  gravity  conditions,  the  Marangoni  force  is  negligible  com¬ 
pared  to  the  buoyant  force.  On  the  other  hand,  under  orbital  conditions 
where  the  maximum  effective  gravity  is  less  than  10-6  g,  the  Marangoni 
force  dominates  for  dT/dz  “  1.82  K/m  (l°F/ft).  Assuming  the  bubble 
motion  is  according  to  Stokes*  law,  the  drag  force  on  the  bubble,  F^  is 


F  =  6pm  R  V 
d  o  0 


and  the  velocity  of  a  bubble  under  the  influence  of  a  surface  tension 
gradient  for  a  zero-gravity  case  can  be  estimated  by  equating  F^  and  Fffl, 
giving 
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Iq-2  (2.6  x  10-6  ■j- n v,' ) .  Assuming  R  «•  0.30  cm,  the  bubble 
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velocity  due  to  a  0.182  K/cm  (10°F/ft)  temperature  gradient  is 
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From  the  analysis  it  can  be  concluded  that  for  purposes  of  heat  transfer 
analysis  in  the  zero-gravity,  zero  contact  angle  case,  the  bubble  can  be 
assumed  to  be  located  at  the  hottest  point  in  the  liquid. 
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3.3.4  Boundary  layer  effects 


Boundary  layer  effects  in  the  thin  layer  of  liquid  near  the  point 
where  the  liquid-vapor  interface  contacts  the  wall  is  another  aspect  of 
void  behavior  that  nust  be  understood  because  of  the  important  part  this 
region  plays  in  solidification.  In  order  to  simplify  the  analysis,  the 
system  considered  for  the  boundary  layer  problem  will  consist  of  a  flat 
plate  with  an  axisymmetric  bubble  contacting  the  plate  from  below.  The 
contact  angle  is  zero  and  gravity  is  1-g. 

To  further  specify  the  problem,  the  bubble  volume  is  set  equal  to 
that  of  a  0.254-cm  radius  sphere.  This  system,  with  bubble  shape  calcu¬ 
lated  using  Eqs.  (3.5)  and  (3.6),  is  shown  in  Fig.  3.11. 

In  the  region  near  the  wall  the  interface  curvature  is  small,  so  a 
relatively  simple  equation  relating  the  pressure  difference  across  the 
interface  to  interface  shape  can  be  developed.  The  derivation  follows 
the  method  given  in  Ref.  13.  The  variation  of  the  free  energy,  F,  of  an 
isothermal  capillary  system  is  given  by 


<SF  —  ( P  -  P„)  <5V„  +  a6A 
V  SL  l  c 
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and  the  area  of  the  deformed  element,  dA,  is 


dA  = 


f(§) 


dA 


where  h  is  the  vertical  displacement  from  the  plane  z 
values  of  h,  the  above  equation  can  be  approximated  by 


=  0.  For  small 
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the  change  in  the  liquid  volume  produced  by  a  small  displacement  of  the 
surface  is 


6V  «  2 it  /  r  <5h  dr  . 
i  J  o 


(3.10) 


Substituting  Eqs.  (3.9)  and  (3.10)  into  Eq.  (3.8),  the  change  in  the 
free  energy  is 
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For  an  equilibrium  state  $F  ■  0.  Since  5h  is  arbitrary,  it  follows  that 
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(3.11) 


The  approximation  implied  by  Eq,  (3.11)  can  be  determined  by  comparing 
Eq.  (3.11)  to  Eq.  (3.3).  Equation  3.3  can  be  rewritten  as 


P  -  P. 

v  l 


__  A  dz  1  dr 

r  dll  dL 


Near  the  point  of  contact  with  the  horizontal  plate  dr  +  dL.  Making 
this  substitution  in  the  above  equation  gives 


p  -  P. 
v  i 


which  when  rearranged  is  identical  to  Eq.  (3.11).  Thus,  the  approxima¬ 
tion  made  in  Eq.  (3.11)  is  that  the  interface  length  dL  a  dr.  If  it  is 
further  assumed  that  the  difference  in  hydrostatic  head  can  be  neglected 
across  the  boundary  layer,  the  pressure  difference  is  constant  and 
Eq.  (3.11)  can  be  solved  analytically.  The  resulting  equation  for  h  as 
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a  function  of  r  is 
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(3.12) 
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V  A 

where  c  =  - - - ,  a  constant,  and  rc  is  the  radius  at  which  the  liquid- 

vapor  interface  contacts  the  wall.  (—  =  0  at  r  =  r  was  used  as  a 

dr  c 

boundary  condition.)  Figure  3.12  shows  the  small  displacement  approxi¬ 
mation  solution,  Eq.  (3.12),  compared  to  the  exact  solution  (also  shown 
in  Fig.  3.11)  for  the  example  case.  The  results  plotted  in  Fig.  3.12 
indicate  that  the  small  displacement  approximation  agrees  reasonably 
well  with  the  exact  solution  for  h  <  0.0254  cm. 

In  order  to  determine  the  effect  of  solidification  in  the  boundary 
region,  the  momentum  integral  equation  can  be  applied  over  the  boundary 
layer  to  derive  a  differential  equation  for  the  boundary  layer  thick¬ 
ness,  h  (r).  The  momentum  integral  equation  for  the  present  case  can  be 
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Fig.  3.12.  Comparison  of  small  displacement  solution  with  exact 
solution. 


derived  by  applying  the  conservation  of  momentum  principle  to  a  control 
volume  as  shown  in  Fig.  3.13  (Ref.  14).  Since  the  flow  is  axisymmetric, 
the  control  volume  forms  an  annular  ring  about  the  axis  of  symmetry. 
Total  radial  momentum  entering  and  exiting  the  control  volume  is  found 
by  integrating  over  the  entire  boundary  layer  thickness.  The  change  in 
momentum  is  equated  to  the  sum  of  forces  on  the  control  volume,  includ¬ 
ing  viscous  stresses  and  pressure.  The  resulting  equation  is 


2ir r  /  pu2dz  .  —  2irr  pu2dz 

o  r  —  &r  1  o  r 


2lT r  3h<P|  r-F  r_Ar)  +  Ar  (Th~  V3 


where  u  is  the  radial  velocity  component,  t  is  the  shear  stress  at  the 

o 

wall  and  is  the  shear  stress  at  the  free  surface.  After  rearranging 
and  letting  Ar  0,  t\e  final  form  of  the  momentum  integral  equation  is 
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Fig.  3.13.  Control  volume  for  momentum  integral  analysis. 
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In  order  to  solve  Eq.  (3.13)  a  velocity  profile  must  be  assumed.  For 
creeping  flow  along  a  plate  the  velocity  profile  is  parabolic,  so  the 
velocity  profile  is  assumed  to  be  of  the  form 


u-A+B—  +  C 
h 


where  A,  B,  and  C  are  constants.  Boundary  conditions  on  u  are 


u  =  0  at  z  “  0 


i,  at  2  *  h 
n 


Applying  the  boundary  conditions,  the  velocity  profile  can  be  rewritten 


u  =  B 


©  •«(--)  («) 


The  value  of  the  constant  B  is  determined  by  mass  conservation.  Solid¬ 
ification  is  modeled  by  a  constant  mass  flux  per  unit  area  through  the 
wall,  m.  Assuming  the  boundary  layer  extends  to  the  centerline,  as 
shown  in  Fig.  3.13,  the  requirement  for  conservation  of  mass  can  be 
written 

*  9  rh 

m  •  irr^  =  2irr  J  pudz 
o 

Inserting  the  expression  for  the  velocity  profile  into  the  above  equa¬ 
tion,  the  value  of  the  constant  is  found  to  be 

xuh  o  • 

B  + 

2p  2  ph 

and  the  velocity  profile  is  therefore 
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The  shear  stress  at  the  wall  can  be  found  by  differentiating  Eq.  3.14, 
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(3.15) 


Setting  Py  =  0  in  Eq.  (3.11)  the  pressure  in  the  liquid  is  given  by 


(3.16) 


Eqs.  (3.14),  (3.15),  and  (3.16)  can  now  be  inserted  into  Eq.  (3.13)  and 
after  much  algebra,  the  resulting  differential  equation  for  h  is 
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Equation  (3.17)  was  solved  by  the  Runge-Kutta  method  using  repre- 

sentative  values  of  m  and  x,  .  It  was  assumed  that  xL  resulted  from  a 

h  h 

surface  tension  gradient.  Boundary  conditions  were  provided  by  matching 
with  the  static  solution  given  by  Eq.  (3.12)  at  h  =  0.0254  cm.  It  was 
found  that  the  deviation  from  the  static  solution  was  negligible  for 
h  >  2.54  x  10-l+  cm.  This  result  indicates  that  the  static  interface 
shape  given  by  Eq.  (3.5)  and  Eq.  (3.6)  is  approximately  correct  in  the 
boundary  region  during  solidification. 

In  summary,  the  analysis  indicates  that  in  microgravity  the  void 
will  be  located  at  the  hottest  point  in  the  liquid,  and  thus  In  the  case 
of  uniform  cooling  a  concentric  void  will  form.  An  analysis  based  on 
capillary  hydrostatics  appears  to  model  void  behavior  in  an  axisvmmetric 
solidification  problem  with  sufficient  accuracy.  However,  it  should  be 
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noted  that  in  nany  instances  liquid  lithium  hydride  has  been  observed  to 
"crawl"  up  the  walls  of  containers  in  a  manner  that  seems  inconsistent 
with  what  the  hydrostatic  analysis  would  predict.  The  theory  to  explain 
such  behavior  is  not  available  at  present.  Hopefully,  an  understanding 
of  this  behavior  will  be  gained  from  the  experiments  to  be  conducted  as 
part  of  this  program. 

3.4  Phase-Change  Induced  Stress 

3.4.1  Spherical  shell  elastic  stress  analysis 

As  indicated  in  Sect.  3.2,  the  void  is  expected  to  be  located  in 
the  sphere  center  following  cooldown  in  a  micro-gravity  space  environ¬ 
ment.  Then,  during  heat  up,  the  liquid  lithium  hydride  expands  against 
the  shell  and  solid  lithium  hydride  as  shown  in  Fig.  3.14a.  Figure 
3,14b  shows  the  results  of  a  thermal  analysis  performed  for  a  3.81-cn 
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Fig.  3.14a.  Heat-up  scenario  of  spherically  encapsulated  lithium 
hydride  in  micro-gravity. 
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Fig.  3.14b.  Spherically  encapsulated  lithium  hydride  thermal 
profile  during  heat-up,  Rq  =  1.9  cm. 
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diam  sphere  with  an  initial  temperature  of  700  K  suddenly  exposed  to 
liquid  netal  at  1100  K.  The  one-dimensional  "crush"  model  described  in 
Sect.  3.2  was  used  to  calculate  the  shell  temperature,  average  tempera¬ 
ture  of  solid  lithium  hydride  and  minimum  temperature  of  solid  lithium 
hydride  as  a  function  of  time.  It  can  be  seen  that  while  the  shell 
reaches  the  liquid  metal  temperature  very  rapidly,  much  of  the  solid  is 
at  a  temperature  significantly  below  the  melt  point  when  the  lithium 
hydride  adjacent  to  the  shell  starts  to  melt.  Since  the  shell  has  a 
limited  capacity  to  expand  elastically,  there  are  three  possible  heatup 
scenarios : 

1.  rupture  (plastically  deform)  the  encapsulating  shell, 

2.  "break."  (deform  or  collapse)  the  lithium  hydride  shell,  or 

3.  liquid  leaks  through  the  solid  lithium  hydride  (via  cracks)  to  the 
void. 

If  it  is  assumed  that  the  void  forms  in  the  sphere  center,  and  the 
lithium  hydride  contains  no  cracks,  then  a  simple  elastic  analysis  can 
be  performed  to  determine,  in  an  approximate  manner,  the  required  shell 
thickness  to  prevent  shell  rupture.  Oliver15  used  elastic  analysis  to 
examine  the  potential  use  of  a  silicon  carbide  shell  for  lithium  hydride 
containment . 

The  required  shell  thickness  for  shells  constructed  of  304L  stain¬ 
less  steel  and  molybdenum  were  determined  using  an  elastic  analysis.  In 
this  analysis  the  following  assumptions  are  made: 

1.  the  void  forms  in  the  sphere  center, 

2.  the  solid  lithium  hydride  forms  a  thick  shell  with  no  cracks, 

3.  the  containment  shell  is  a  thin  shell,  and 

4.  thermal  stresses  other  than  those  resulting  from  phase-change  are 
neglected. 

In  addition,  based  on  the  thermal  profiles  in  Fig.  3.14b,  the  con¬ 
tainment  shell  is  assumed  to  be  isothermal  at  the  peak  operating  tem¬ 
perature,  and  the  solid  lithium  hydride  is  assumed  to  be  isothermal  at 
the  minimum  storage  temperature.  These  assumptions  lead  to  a  conserva¬ 
tively  low  shell  tensile  strength  and  a  conservatively  high  lithium 
hydride  compressive  strength. 


A  force  balance  on  a  tnin-wall  sphere  with  internal  pressure,  P, 
and  no  external  pressure  yields:*6 


2t  °s/r  * 


(3. 18) 


where 

t  “  shell  thickness, 
r  »  spherical  shell  radius,  and 
Os  -  shell  tensile  strength. 

Timoshenko  and  Goodier17  have  developed  the  stress  equations  for  a 
thick-wall  spherical  shell.  The  maximum  principal  stress  occurs  in  the 
tangential  direction  at  the  inner  radius,  and  is  compressive: 


(3.19) 


(r5  -  a3) 


where 


a  ”  maximum  lithium  hydride  compressive  stress, 

LiH 

a  »  inner  radius  of  the  lithium  hydride  shell,  and 

r  «  outer  radius  of  the  lithium  hydride  shell  (encapsulating 

shell  radius). 

As  shown  by  Oliver,15  if  the  sphere  is  assumed  to  be  completely  full  at 

the  hottest  liquid  condition  and  a  central  void  forms  on  cooldown,  then: 


(r/a)3  -  (1  -  P*/Ps>'1 


(3.20) 


where 

P . ,  p  are  the  liquid  and  solid  lithium  hydride  densities  at  the 
it  s 

hottest  temperature  and  minimum  storage  temperature, 
respectively. 

Using  Eqs.  3,18  through  3.20,  it  follows  that  the  minimum  required 
encapsulating  shell  thickness  to  prevent  rupture  is  given  by: 


r/3  '■WV'W 


(3.21) 


Pressed  lithium  hydride  ultimate  compression  strength  data  are  given  in 
Fig.  3.15  (left  vertical  axis).  As  shown  there  is  a  strong  temperature 
dependence.18*13  Table  3.3  presents  lithium  hydride  compressive 
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Fig.  3.15.  Pressed  and  estimated  cast  lithium  hydride  ultimate 
compression  strength  (Refs.  18  and  19). 


Table  3.3.  Compressive  strength  of  lithium  hydride 
at  room  temperature  (Ref.  20) 


Strength 

-  Comments 
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30.13 
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± 

6.31 

24,300 

± 

915 

(95Z  confidence  Interval 
room-temperature  data) 


Cold  pressed 

Cold  pressed  and  sintered,  3  cycles* 

Cold  pressed  and  sintered,  5  cycles* 

Cold  pressed  and  sintered,  10  cycles* 

‘Blocks  therma1  cycled  from  room  tem¬ 
perature  to  593I>C,  specimens 
machined  from  sound,  uncracked 
portion  of  blocks 


strength  at  room  temperature  for  pressed  and  sintered  material.  ^0  These 
data  indicate  that  sintering  results  in  large  strength  gains  (about  a 
factor  of  1.7  greater  for  the  10  cycle  sample).  The  sintered  data  are 
more  likely  representative  of  cast  material,  thus  the  data  of  Fig.  3.15 


should  probably  be  scaled  upward  by  this  factor.  This  has  been  done  in 
the  right  vertical  axis  of  Fig.  3.15. 

At  — 11 00  K,  the  yield  strength  of  304L  stainless  steel21  and 
molybdenum22  are  taken  as  69  (estimated)  and  255  MPa  (10  and  37  ksi), 
respectively.  These  properties  do  not  reflect  any  effects  of  hydrogen 
or  lithium  on  strength;  these  effects  are  discussed  in  Sect.  3.6. 

A  small  computer  program,  SHELL,  was  written  to  calculate  achiev¬ 
able  energy  storage  densities  (i.e.,  includes  shell  and  lithium  hydride) 
based  on  the  minimum  required  shell  thickness  determined  as  described 
above.  Results  are  presented  in  Table  3.4  for  encapsulated  lithium 
hydride  in  spherical  shells  constructed  of  304L  stainless  steel  and 
molybdenum.  Storage  densities  are  based  on  a  500  s  heatup  period  and 
use  of  naturally  occurring  lithium  and  are  shown  for  a  range  of  minimum 
storage  temperatures.  These  results  show  that  the  required  shell  thick¬ 
ness  is  very  large  at  low  minimum  storage  temperatures,  reflecting  high 
lithium  hydride  compressive  strength  at  low  temperatures.  Obviously, 
with  thick  shells,  only  a  very  low  energy  storage  can  be  achieved.  Even 
for  minimum  storage  temperatures  in  the  500  to  700  K  range,  calculated 
shell  thicknesses  are  much  greater  than  permitted  to  obtain  desired 
energy  storage  densities.  As  the  lithium  hydride  melt  temperature  is 
approached  and  the  lithium  hydride  strength  becomes  small,  more  reason¬ 
able  shell  thicknesses  and  energy  storage  densities  can  be  achieved. 
For  comparison  purposes,  energy  storage  with  lithium  is  shown  in  Table 
3.4  where  it  is  assumed  that  a  shell  with  negligible  mass  can  be  used. 
In  the  operating  range  of  interest  for  the  minimum  storage  temperature 
(500  to  700  K) ,  the  achievable  storage  density  is  still  well  below  that 
desired. 

The  storage  densities  presented  in  Table  3.4  can  be  compared  to 
those  shown  in  Table  2.1  for  Li^H  and  Li7H,  alone,  to  see  the  penalty 
paid  for  the  encapsulating  shell  or  using  lithium  instead  of  lithium 
hydride.  Based  on  elastic  analysis  it  appears  that  a  desired  shell 
thickness  of  ~0.0127  cm  (5  mil)  cannot  be  obtained  in  the  temperature 
range  of  interest  (even  for  a  molybdenum  shell,  a  material  with  substan¬ 
tial  high  temperature  strength)  if  lithium  hydride  forms  a  monolithic 
she  1 1 . 
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However,  it  has  been  observed  in  tests  conducted  for  this  program 
and  others , 23 * 24a  that  cracks  form  in  the  lithium  hydride  during  cool¬ 
down  (see  Sect.  4.1).  Thus,  cracks  that  penetrate  the  solid  lithium 
hydride  from  the  containment  shell  to  the  void,  may  provide  a  path  for 
expanding  liquid  lithium  hydride  to  reach  the  centrally  located  void 
during  heatup.  Hence,  large  hydrostatic  forces  are  avoided,  and  a  thin 
shell  can  successfully  contain  the  lithium  hydride. 

In  the  event  that  predictable,  "well-behaved"  cracks  do  not  form, 
there  are  a  series  of  potential  design  modifications  to  mitigate  the 
stress  problem.  These  are  discussed  in  the  following  section. 


induced  stress 


As  previously  mentioned,  the  potential  for  high  stresses  in  the 
shell  resulting  from  melting  lithium  hydride  when  the  liquid  is  not  in 
contact  with  the  void  is  one  of  the  most  serious  potential  problems  in 
developing  a  durable  energy  storage  capsule.  Design  modifications  that 
have  been  investigated  as  potential  means  of  reducing  phase-change 
induced  stress  are  listed  below: 

1.  providing  internal  fins  made  of  the  shell  material 

2.  making  the  container  from  a  non-wetting  material  or  insulating  part 
of  the  surface,  thus  causing  the  void  to  form  at  the  wall 

3.  making  the  container  flexible 

These  potential  modifications  are  discussed  below. 

Internal  fins 

Two  potential  benefits  result  from  providing  internal  fins  that 
extend  from  the  shell  into  the  solid  lithium  hydride.  First,  they  keep 
the  lithium  hydride  from  forming  a  monolithic  solid,  thus  making  the 
solid  easier  to  break  (the  concept  is  somewhat  analogous  to  removing  ice 
cubes  from  an  ice  tray).  In  addition,  by  providing  a  thermal  conduction 
path  to  the  void,  the  fins  provide  a  potential  channel  for  liquid  to  bv- 
pass  the  solid  into  the  void.  The  effectiveness  of  fins  as  a  heat  con¬ 
duction  path  can  be  estimated  by  a  constant  density  finite  difference 
calculation.  For  the  example  problem  shown  in  Fig.  3.16,  lithium 
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Fig.  3.16.  Spherically  encapsulated  lithium  hydride  with  an 
internal  fin  shown  in  cross-section  in  left  symmetric  half. 


hydride  in  a  1.9-cm  radius  0.0381-cm  (15-mil)  thick  spherical  container 
is  considered.  A  fin  made  of  the  shell  material  extends  completely 
around  the  equator  of  the  shell.  Figure  3.17  shows  the  effect  of  the 
fin  on  the  temperature  distribution  in  the  solid  lithium  hydride  when  a 
sphere  with  initial  temperature  of  700  K  (800°F)  is  suddenly  exposed  to 
convection  from  NaK  at  1100  K  (1520°F).  Temperature  distributions  in 
the  lithium  hydride  adjacent  to  the  fin  and  along  the  centerline  (at  a 
90°  angle  from  the  fin)  after  1  s  of  heating  are  shown  for  shells  made 
of  stainless  steel,  niobium,  and  copper.  (It  should  be  noted  that  these 
temperatures  are  average  values  for  the  finite  difference  control  vol¬ 
umes  adjacent  to  the  fin;  temperatures  at  the  fin  surface  are  higher.) 


The  effect  of  the  fin  is  small  in  the  case  of  stainless  steel  and 


niobium.  Only  copper  has  a  conductivity  sufficiently  higher  than  the 
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Fig.  3.17.  Comparison  of  temperature  distri 
Lithium  hydride  adjacent  to  fin  versus  90°  away  f 
hea  ting. 


lithium  hydride  to  provide  a  conduction  path  tha 
the  mechanical  strength  of  the  lithium  hydride  ni 
mal  d  i  f f us i vi t ies  of  stainless  steel,  niobium  a 
tively,  about  5,  25,  and  120  times  that  of  soli 

can  be  concluded  from  the  above  analysis  that 
fusivity  would  need  to  be  about  two  orders  of 
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that  of  solid  lithium  hydride  in  order  for  fins  to  be  an  effective 
conduction  path.  It  is  not  known  if  potential  shell  materials  with  such 
high  thermal  dif fusivit ies  exist.  However,  even  without  a  significant 
heat  conduction  effect  fins  would  be  effective  if  the  solid  lithium 
hydride  does  not  form  a  strong  bond  with  the  shell  material  and  is  a 
technique  that  merits  further  investigation. 


Void  management 


In  a  micro-gravity  environment,  the  liquid-vapor  interface  is  a 
spherical  segment  with  the  radius  being  determined  bv  the  contact 
angle.  For  a  non-zero  contact  angle  the  void  would  be  most  likelv  to 
form  at  a  wall,  since  a  wall-bound  bubble  has  less  surface  area  per  unit 
volume  than  a  free  bubble.**  Therefore,  it  would  seem  that  if  the 

container  is  made  of  a  material  which  lithium  hydride  does  not  wet  the 

void  would  form  at  the  wall.  It  is  not  known  at  this  time  if  such 

materials  exist  and  efforts  are  underway  to  obtain  data  on  lithium 
hydride  wetting  behavior.  In  any  event,  a  cold  container  surface 

exposed  to  lithium  hydride  vapor  would  tend  to  form  a  thin  surface  layer 
of  frozen  lithium  hydride  as  a  result  of  condensation.  The  void  would 
then  assume  the  zero  contact  angle  orientation.  Therefore,  it  does  not 
seem  likely  that  wetting  behavior  alone  can  be  used  to  control  the  void 
location. 

Another  possible  method  for  influencing  void  location  would  he  to 
reduce  the  cooling  rate  on  part  of  the  container.  This  would  tend  to 
cause  the  void  to  form  at  the  area  of  slowest  cooling  and  might  result 
in  a  thinner  layer  of  solid  being  formed  in  that  area.  A  full  evalua¬ 
tion  of  this  idea  will  require  a  heat  transfer  computer  code  with  an 
integrated  void  mechanics  model.  The  influence  of  insulating  a  small 
section  of  the  container  surface  on  the  rate  of  solidification  can  be 
demonstrated  by  a  constant  density  finite  difference  calculation.  For 
the  example  problem  shown  in  Fig.  3.1b,  the  surface  of  a  1.9-cm  radius 
sphere  containing  liquid  lithium  hydride  is  exposed  to  liquid  metal 
convection  except  for  an  area  extending  13°  from  the  axis  ot  symmetry  on 
one  end,  which  is  insulated.  The  temperature  of  the  heat  transfer  fluid 
decreases  gradually  with  time,  si  mu  1  iting  rejection  of  heat  to  space 
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Fig.  3.18.  Effect  of  an  insulating  surface  on  the  solidification 
of  spherically  encapsulated  lithium  hydride. 


through  a  radiator.  As  shown  in  the  figure,  the  thickness  of  the  frozen 
layer  under  the  insulated  area  is  much  less  than  under  the  uninsulated 
surface  after  cooling  for  800  s,  and  this  effect  should  be  Increased 
when  the  influence  of  the  void  is  taken  into  consideration.  The  dashed 
line  in  Fig.  3.18  shows  the  volume  a  concentric  void  would  occupv  when 
the  lithium  hydride  is  completely  solid. 

Unfortunately,  it  will  not  be  possible  to  evaluate  void  management 
schemes  experimentally  under  normal  gravity  conditions.  Experiments  in 
orbit  will  be  needed  to  confirm  the  predicted  performance. 


Flexible  containers 


Phase-change  induced  stress  in  the  container  walls  can  be  elimi¬ 
nated  if  the  container  can  be  designed  to  accommodate  the  excess  liquid 
volume  produced  by  melting  via  expansion.  Estimates  of  the  maximum 
volume  expansion  that  axially  symmetric  containers  can  undergo  without 
damage  can  be  made  by  considering  the  container  to  be  a  thin  elastic 
membrane.  Equations  for  the  deformations  of  shells  of  revolution 
subjected  to  uniform  internal  pressure  are  given  in  Ref.  25.  These 
equations  were  solved  for  a  sphere  and  for  oblate  spheroids  (ellipse 
rotated  about  its  minor  axis)  of  various  eccentricities.  The  flattest 
spheroidal  shell  that  can  be  subjected  to  internal  pressure  without 
buckling  has  a  ratio  of  semiaxes  =  /T.  It  was  found  that  for  oblate 

spheroids  with  —  <  /2  a  sphere  (—  =  1)  has  the  greatest  volume  expansion 
b  b 

for  a  given  wall  thickness  and  maximum  allowable  stress  level.  In  anv 
case,  the  maximum  volume  expansion  that  can  he  achieved  with  a  practical 
shell  design  is  about  1%,  which  is  far  less  than  needed  to  accommodate 
the  liquid  volume  expansion. 

On  the  other  hand,  bellows  containers  can  acconnodate  rather  large 
volumetric  expansions,  on  the  order  of  20  to  30?,  without  incurring 
damage.  Figure  3. 19  shows  a  photograph  of  a  304L  stainless  steel  seam¬ 
less  bellows  with  dimensions  in  the  range  of  interest:  length  of 
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Fig.  3.19.  304L  stainless  steel  bellows  container. 
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As  shown  in  Fig.  3.20,  the  required  hydrogen  overpressure  to 
prevent  lithium  hydride  dissociation  for  high  mole  fraction  lithium 
hydride  systems  increases  dramatically  with  increasing  lithium  hydride 
mole  fraction.  For  a  pure  lithium  hydride  system  (LiH  mole  fraction  of 
100%),  the  equilibrium  hydrogen  pressure  would  theoretically  be 
infinite.26  Thus,  for  a  practical  system  it  will  be  necessary  that  the 
encapsulat-  d  lithium  hydride  contain  excess  lithium  to  permit  a  reason¬ 
able  hydrogen  overpressure. 

Since  lithium  hydride  tends  to  dissociate,  free  hydrogen  is  present 
in  the  encapsulating  container  and  can  diffuse  through  the  walls.  The 
loss  of  hydrogen  from  the  container  results  in  a  loss  of  energy  storage 
density.  As  described  in  Ref.  27,  hydrogen  diffusion  (loss)  can  be 
calculated  using  the  following  equation: 


H2  loss 


S  t  A  P1/2 


(3.22) 


where 

S  =  hydrogen  permeability  of  the  metal, 
t  =  time, 

A  =  surface  area  of  the  encapsulating  metal, 

P  =  hydrogen  pressure,  and 
d  =  shell  wall  thickness. 

Permeabilities  of  stainless  steel  and  molybdenum,  the  materials  of  in¬ 
terest  for  encapsulating  lithium  hydride,  are  shown  in  Table  3.5 
(Ref.  27)  along  with  permeabilities  of  other  materials.26  As  shown,  the 
permeability  of  clean  stainless  steel  is  quite  high  whereas  the  permea¬ 
bilities  of  oxidized  stainless  steel  and  molybdenum  are  nuch  lower. 

Hydrogen  loss  from  lithium  hydride  encapsulated  in  molybdenum  and 
clean  stainless  steel  has  been  calculated  for  a  25  cycle  lifetime.26 
The  calculations  were  performed  for  a  sphere  with  a  diameter  of  1.5  in. 
containing  approximately  15.5  grams  of  LiH  (initial  mole  fraction 
99.16%).  In  addition,  calculations  assumed  the  use  of  naturally  occur¬ 
ring  isotopic  hydride,  a  thermal  storage  temperature  range  of  700  to 
1100  K,  and  a  pulse  duration  and  orbit  time  of  500  and  5880  s,  respec¬ 
tively.  Phenomenologically,  it  was  assumed  that  hydrogen  diffusion 
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Table  3 

.5.  Hydrogen  permeabilities 

of 

selected  materials27 

(Pressure  =  760  mmHg) 

Material 

„  cc( STP)mm 

Permeability  - r-r- - ; 

h-atm1'  '-cm^ 

SS  304  -  clean 

0.4  5a 

SS  304  -  oxidized 

0 . 0 1  2a 

Molybdenum 

0.07a 

SiC 

4.3  x  10'6* 

Tungsten 

0.0007* 

Nickel 

1.18* 

Iron 

1.47* 

Temperature  near  1100  K. 
''Temperature  near  1000  K. 


occurred,  as  calculated  by  Eq.  3.22,  when  liquid  lithium  hydride  or 
hydrogen  gas  contacted  the  container  wall,  and  that  on  cooldown  lithium 
hydride  freezes  first  and  uniformly  on  the  shell  surface  preventing 
further  hydrogen  diffusion  (i.e.,  the  lithium  hydride  freezes  over  the 
entire  shell  surface  which  ultimately  results  in  a  central  void).  The 
pressure  used  in  Eq.  3.22  was  allowed  to  vary  from  cycle  to  cycle  and 
was  assumed  to  be  the  hydrogen  equilibrium  pressure  corresponding  to  the 
lithium  hydride  composition  at  the  beginning  of  each  cycle. 

Based  on  the  assumptions  described  above,  Table  3.6  presents  the 
calculated  loss  in  energy  storage  density  for  molybdenum  and  clean 
stainless  steel  shells  0.00254,  0.0127,  and  0.0381  cm  (1,  5,  and  10  mils) 
Chick.  These  results  indicate  that  hydrogen  loss  through  a  0,0127  cm 
(5  mil)  clean  stainless  shell  is  only  2.9%;  however,  the  loss  from  a 
0.00254  cm  (1  mil)  clean  stainless  steel  shell  is  unacceptably  high  at 
11%.  On  the  other  hand,  molybdenum  shells  as  thin  as  0.00254  cm  (1  mil) 
are  acceptable  with  less  than  3%  loss  in  storage  density.  Since  the 
hydrogen  permeability  of  oxidized  stainless  steel  is  significantly 
better  than  molybdenum  values,  a  0.00254  cm  (1  mil)  oxided  stainless 
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Table  3.6.  Effect  of  hydrogen  loss  on 
energy  storage  density 


Shell  thickness 
(cm)  (mil) 

Material 

Energy  density 
(MJ/kg) 

1st  cycle  25th  cycle 

Energy  density 
loss  (Z) 

0.00254 

1 

SS-clean 

5.468 

4.880 

11. 

0.0127 

5 

SS-clean 

4.493 

4.362 

2.9 

0.0381 

15 

SS-c lean 

3.112 

3.068 

1.4 

0.00254 

1 

Mo  lybdenum 

5.379 

5.250 

2.4 

0.0127 

5 

Molybdenum 

4.206 

4.173 

0.78 

0.0381 

15 

Molybdenum 

2.725 

2.715 

0.37 

steel  shell  would  certainly  be  acceptable.  The  nominal  stainless  steel 
oxide  coating  would  not  be  stable  in  this  application.  Thus,  use  of  a 
very  thin  stainless  steel  shell  (~1  mil)  would  require  identification  of 
a  stable,  low  hydrogen  permeability  coating  (possibly  some  type  of 
refactory  ceramic)  to  apply  to  the  outside  shell  surface.28 

3.6  Material  Considerations 

Successful  encapsulation  of  lithium  hydride  requires  that  the  shell 
material  be  compatible  with  lithium  hydride,  lithium  and  hydrogen  since 
all  three  materials  will  be  present  in  the  capsule.  Compatibility  is 
required  over  the  temperature  range  ~700  to  1100  K,  where  most  of  the 
exposure  is  at  ~700  K.  As  mentioned  previously,  since  the  system  is 
expected  to  be  cycled  25  times  over  its  lifetime,  relatively  short  life 
is  required  at  high  temperatures  (I.e.,  above  the  lithium  hydride  melt 
temperature ) . 

In  addition  to  being  compatible  with  lithium  hydride,  lithium,  and 
hydrogen,  the  shell  must  be  resistant  to  thermal  shock,  since  the  heatup 
rate  is  severe.  Thus,  the  shell  material  must  be  reasonably  ductile 
with  sufficient  strength  to  contain  lithium  hydride  in  a  thin  wall  shell 


in  light  of  thermal  str  sses,  phase-change  induced  stress  on  heatup,  and 
solidification  shrinkage  forces  generated  during  cooldown.  In  addition, 
material  costs  must  be  considered  as  well  as  the  relative  ease  of 
fabricating  capsules. 

The  results  of  lithium  hydride  compatibility  studies  have  been  re¬ 
ported  in  several  sources. 30-35  Figure  3.22  presents  results  of  stress 
corrosion  testing  of  a  group  of  metals  in  contact  with  lithium  hydride 
at  866  K.30  as  shown,  only  430  stainless  steel  survived  the  full  1000  h 
without  failure;  however,  the  material  has  relatively  low  strength  at 
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Fig.  3.22.  Results  of  stress  corrosion  testing  of  various  metals 
in  contact  with  lithium  hydride  at  866  K. 


Fig.  3.22  source: 


Hamill,  C.  W.,  "Stress  Corrosion  Testing  of  Various 
Alloys,"  Union  Carbide  Nuclear  Company,  in  "Lithium 
Hydride  as  a  Thermal  Energy  Storage  Material," 
Woods,  Frank  L. ,  Wannemacher,  Marvin  P. ,  and  Houck, 
Oscar  0.,  editors,  Wright-Patterson  Air  Force  Base, 
ASD-TR-6 1-427 ,  May  1962,  Fig.  53,  p.  79. 
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high  temperature.  Uastelloy  B,  Inconel,  Inconel  (exposed  to  lithium 
hydride  with  1%  F  as  LiF),  Inconel  600  and  316  stainless  steel  survived 
for  500  h  or  more  at  866  K.  In  tests  at  644  K,  316  and  304  (with 
1%  boron  by  weight)  stainless  steels  survived  the  full  1000  h.  Stain¬ 
less  steel  316  also  survived  the  full  1000  h  In  tests  at  811  K.  It  is 
reported  that  at  811  K,  301  stainless  steel  performed  the  best. 

In  another  study  of  lithium  hydride  corrosion,31  molybdenum, 
niobium,  Hastelloy  C,  and  Hastelloy  N  samples  were  tested  at  977  K  in 
Hastelloy  capsules.  Molybdenum  was  found  to  be  the  best  material,  while 
niobium  disintegrated  during  testing.  It  was  reported  that  the  niobium 
disintegration  was  thought  to  be  due  to  the  bimetallic  effect  rather 
than  corrosion. 

Royer32  has  indicated  that  pure  iron  (ARMCO)  and  low  carbon,  300 
series  stainless  steel  are  suitable  lithium  hydride  containment  mate¬ 
rials  below  1023  K,  and  that  for  higher  temperatures  tungsten  and 
molybdenum  might  be  used. 

Smith  and  Miser240  have  compiled  a  data  base  of  lithium  hydride 
(and  lithium)  cotnpa t i bi  1  i ty  and  corrosion  information.  Based  on  studies 
by  Messer,33  it  is  reported  that 

"molten  lithium  hydride  of  high  purity  does  not  appear  to 

attack  Armco  iron,  312  or  347  stainless  steels  at  tempera¬ 
tures  below  720°C." 

Messer  does  not  recommend  mild  steel  or  silica  for  molten  lithium 
hydride,  even  though  others  have  reported  favorable  experience  with 
lit  hi un. 

Lithium  hydride  corrosion  studies  have  been  conducted  at  General 
Electric  by  Welch34*35  where  austenitic  stainless  steels,  and  the 
19-9DL  alloy  (Universial  Cyclops  alloy  austenitic  stainless  steel  type) 
were  investigated.  It  was  found  that  no  significant  changes  in  ultimate 
tensile  strength  and  0.2%  yield  strength  occurred  in  301,  304,  316, 

316F.LC,  317,  318,  321,  and  347  stainless  steels  after  exposure  to 

lithium  hydride  at  992  K  for  periods  up  to  100  h.  Results  for  the 
19-90L  alloy  were  very  similar.  Yield  strength  data  are  presented  in 
Fig.  3.23  for  these  metals  for  no  exposure,  and  65  and  100  h  exposure. 
In  general,  the  data  (Tables  3.7  and  3.8  for  common  stainless  steels’) 
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Table  3.8.  727  K  percent  elongation  of 

welded  stainless  steels  after  various 
exposures  to  molten  lithium 
hydride  at  992±14  K 


Type 

specimen 

Ba 

BRfe 

CG 

CR1^ 

Type 

Hrs 

exp 

301 

0 

40.1 

33.1 

32.9 

35.6 

65 

30.8 

14.7 

20.0 

25.2 

100 

23.7 

13.0 

22.8 

21.5 

304 

0 

34.0 

25.8 

29.1 

30.0 

65 

29.2 

21.0 

22.7 

24.7 

100 

27.5 

21.2 

22.8 

24.2 

316 

0 

32.3 

19.0 

28.8 

28.8 

65 

28.5 

18.0 

24.0 

23.5 

100 

19.5 

18.0 

25.3 

20.0 

317 

0 

29.8 

22.5 

22.8 

25.2 

65 

25.3 

19.5 

16.8 

14.5 

100 

22.5 

25.0 

18.5 

18.0 

321 

0 

28.3 

27.0 

25.0 

28.9 

65 

20.3 

19.2 

16.7 

16.7 

100 

22.3 

21.0 

17.7 

18.0 

347 

0 

- 

- 

- 

- 

65 

26.7 

20.0 

17.8 

18.7 

100 

22.7 

18.5 

18.0 

19.5 

aWeld  transverse  to  length  of 
standard  strip  specimen;  0.037-inch-thick 
sheet. 


Weld  transverse  to  length  of 
standard  strip  specimen,  ground  flush; 
0.050-inch-thick  sheet. 

Q 

Weld  longitudinal  to  length  of 
standard  strip  specimen;  0.037-inch-thick 
sheet . 

^Weld  longitudinal  to  length  of 
standard  strip  specimen,  ground  flush; 
0.050-inch-thick  sheet. 

Source:  Welch,  F.  H. ,  "Properties  of 

Lithium  Hydride  V:  Corrosion  of 
Austenitic  Stainless  Steels  in  Molten 
LiH,"  General  Electric  Company  APEX-673, 
August  1961,  Appendix  3,  p.  68. 
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Selected  stainless  steels  and  19-9DL  alloy  were  also  exposed  to 
heat  and  hydrogen  determine  the  effect  of  hydrogen  on  mechanical 
properties.  For  these  metals,  the  hydrogen  exposure  caused  mechanical 
property  changes  of  the  same  magnitude  to  occur.  Thus,  Welch  concluded 
that  exposure  to  molten  lithium  hydride  does  not  have  a  significant 
effect  on  the  ultimate  tensile  strength,  and  percent  elongation  of  the 
metals  tested  over  the  tine  periods  and  temperatures  examined. 

For  the  300  series  stainless  steels  Welch  reports, 34a 

"Metal lographic  examination  of  unwelded  and  welded  speci¬ 
mens  of  the  eight  alloys  showed  that  exposure  to  molten 
lithium  hydride  up  to  100  h  at  1325°F  or  65  h  at  1425°F  did 
not  produce  any  significant  attack  (0.001  in.  or  less)." 

At  1047  K.  (1425°F),  corrosive  attack  became  significant  for  100-h 
exposure  to  molten  lithium  hydride  for  304,  317,  and  321  stainless 
steels  (0.005-0.01  cm),  and  showed  a  slight  increase  in  316,  316ELC,  and 
347  stainless  steel.  301  stainless  steel  was  essentially  unaffected  at 
1047  K  (1425°F).  For  the  19-9DL  alloy,  Welch36  found  the  lithium 
hydride  attack  to  be  slight  with  no  definite  indication  of  intergranular 
attack.  Stainless  steel  and  the  19-9DL  alloy  welds  showed  a  slightly 
greater  corrosive  attack  than  the  bare  metal. 

Welch,36  based  on  his  studies  of  the  effects  of  molten  lithium 
hydride  on  300  series  stainless  steels  and  19-9DL  allov,3L,»36  concluded 
that  the  19-9DL  alloy  was  definitely  superior  to  the  others  tested. 

Several  studies  have  examined  the  compatibility  of  various 
materials  with  lithium.36-33  Figure  3.24  presents  a  summary  of  resis¬ 


tance  of  various  materials  to  lithium. 


As  indicated,  pure  iron, 


ferritic-chromium  stainless  steel,  columbium,  tantalum,  and  molybdenum 
show  good  resistance  to  lithium  up  to  at  least  800°C  (application  peak 
design  temperature  is  1100  K  or  B27°C).  However,  since  the  system  only 
operates  for  relatively  short  times  at  peak  temperatures,  it  is  probably 
reasonable  to  consider  austenitic  Cr-Ni  stainless  steel  as  being  accept¬ 
able. 

Reference  37  summarizes  the  results  of  numerous  corrosion  studies 
of  refractory  metals  and  alloys,  among  other  materials.  Studies 
indicate  that  Mo-30  W,  Mo-0.5  Ti  ,  and  Mo-0.5  Ti-Q.OH  Zr  resist  static 
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static  lithium  at  1089  K  for  100  h  is  provided  in  Fig.  3.25  (Ref.  38). 
Of  those  studied,  TiC,  ZrC,  and  Cr3C2  show  good  lithium  resistance. 

The  effect  of  hydrogen  on  the  mechanical  properties  of  many  of  the 
refractory  metals  prevents  their  use.  Niobium  suffers  severe  hydrogen 
embrittlement  when  exposed  to  hydrogen  gas  in  the  500  to  600°C  (773  to 
873  K)  range.33  The  material  loses  its  ductility,  and  about  half  of  its 
strength.  Hydrogen  also  adversely  affects  the  mechanical  properties  of 
tantalum,43  t  itanium, 4  *  > 4:3  and  zirconium.43  On  the  other  hand,  rhenium 
has  been  reported44  to  behave  well  in  hydrogen  at  elevated  temperatures 
while  tungsten  is  reported43  as  not  reacting  with  hydrogen. 

Westphal  and  Worzala40  describe  two  forms  of  hydrogen  damage  in 
steels:  hydrogen  attack  and  hydrogen  embrittlement.  Hydrogen  attack  is 
distinguished  from  hydrogen  embrittlement  in  that  in  the  former,  a 
methane  reaction  occurs,  but  not  in  the  latter.  In  addition,  hydrogen 
embrittlement  is  reversible  and  occurs  below  473  K,  while  hydrogen 
attack  is  usually  considered  irreversible  and  occurs  only  above  473  K. 
In  a  summary  paper  by  Morris47  it  is  reported  that  for  high  temperature 
conditions  (greater  than  866  K),  steels  containing  greater  than  12 7 
chromium,  and  the  austenitic  stainless  steels  are  resistant  to  hydrogen 
attack  in  all  known  applications. 

As  described  previously,  Welch34*33  studied  the  effects  of  heated 
hydrogen  on  the  mechanical  properties  of  selected  stainless  steels  and 
19-9DL  alloy.  For  316  and  318  stainless  steels,  and  19-9DL  alloy,  Welch 
concluded  that  the  losses  in  ultimate  tensile  strength,  yield  strength, 
and  percent  elongation  were  due  primarily  to  hydrogen.  Welch,  in 
arriving  at  this  conclusion,  was  aware  of  reported  stability  of  high 
chromium  content  austenitic  steels  to  hydrogen.33  It  is  interesting  to 
note  that  301  and  316  stainless  steels  (which,  as  indicated  previously, 
suffered  more  significant  reductions  in  percent  elongation  following 
exposure  to  molten  lithium  hydride)  have  less  chromium  content  Mf-l^'M 
than  the  other  common  steels  tested  (SS  304  and  317  .and  a )  ]  <>v , 
18—207;  and  SS  321  and  347,  1 7—197). 


Based  on  the  compatibility  of  shell  materials  with  lithium  hvdride, 
lithium,  and  hydrogen  in  the  temperature  range  of  Interest,  manv  poten¬ 
tial  materials  considered  here  can  he  eliminated  for  possible  use.  Most 
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Fig.  3.25.  Corrosion  resistance  of  various  ceramics  in  static 
lithium  for  100  h  at  1089  K, 


Fig.  3.25  source: 


Cook,  W.  H.  ,  "Corrosion  Resistance  of  Various 
Ceramics  and  Cermets  to  Liquid  Metals,"  Oak  Ridge 
National  Laboratory,  ORNL-2391,  May  31,  I960,  Fig. 
10,  p.  14. 


of  the  non- metals  can  be  eliminated  due  to  incompatibility  with  lithium 
hydride  and/or  lithium.  Exceptions  include  ZrC,  TiC,  and  Cr3C2-  With 
the  exception  of  the  refractory  metals,  nonferrous  metals  are  not  com¬ 
patible  with  lithium.  And,  of  the  refractory  metals,  niobium,  zir¬ 
conium,  titanium,  and  tantalum  suffer  from  hydrogen  embrittlement.  Of 
the  ferrous  metals,  pure  iron,  ferritic-chromium  stainless  steels,  and 
austenitic  Cr-Ni  stainless  steels  offer  either  good  or  limited  lithium 
combatibility  at  temperatures  at  least  as  high  as  800°C. 

In  addition  to  material  compatibility,  other  material  considera¬ 
tions  applicable  to  the  encapsulating  application  include: 

1.  strength, 

2.  density, 

3.  ductility, 

4.  thermal  shock  resistance,  and 

5.  material  and  fabrication  cost. 

Table  3.9  presents  the  yield  strength,  density,  and  ductility  (total 
elongation)  of  materials  (data  from  Refs.  21,  22,  48—53)  which,  based 
on  material  compatibility,  could  be  considered  for  use.  As  indicated, 
pure  iron  and  ferritic-Cr  stainless  steels  have  low  strength  at  high 
temperature.  Unlike  SiC,  TiC  (and  also  ZrC  and  CrjC2)  do  not  have  good 
thermal  shock  resistance. 52  Tungsten  (and  rhenium)  have  high  densi¬ 
ties.  304L  stainless  steel  has  moderate  high  temperature  strength  and 
density  and  good  ductility.  Similarly,  molybdenum  has  moderate  density, 
and  good  ductility  (however,  molybdenum  has  a  high  brittle-to-ductile 
transition  temperature,  ~100°C),  but  also  possess  high  strength  at 
elevated  temperatures. 

It  is  felt  that  constructing  and  sealing  a  shell  using  one  of  the 
carbide  ceramics  would  be  quite  difficult.^4  Fabricating  shells  of 
tungsten  or  rhenium  would  be  more  difficult  than  molybdenum,  which  is  a 
lower  cost  refactory  metal  which  can  be  spun  or  hydrof orned . 55 
Obviously,  the  material  and  fabrication  costs  of  stainless  steels  are 
much  less. 

Based  on  the  data  presented,  it  is  concluded  that  molybdenum  or  one 
of  the  austenitic  stainless  steels  would  provide  the  best  possible 
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Table  3.9.  Selected  mechanical  properties 
of  potential  shell  materials 


Material 

Yield 
st  rength^ 
(MPa) 

Density^ 
(g/cm3 ) 

Total  ^ 

elongat  iorr 
(%) 

Pure  iron48 

-39 

7.87 

26.5-53.5 

aSS-44621.49a.50 

-29 

7.6 

25.h 

^SS-304L21»  49b>  51 

-69 

8.03 

35.-60. 

cTiC52 

296^ 

4.92 

brittle 

Molybdenum^ 

255 

10.2 

24.-26. 

Tungsten22 >  8  8 

-110 

19.3 

r* 

O 

vXD 

1 

• 

o 

^Representative  of  ferritic-Cr  stainless  steels. 

°Representat ive  of  austenitic  stainless  steels. 

^Representative  of  ZrC  and  Cr3C2« 

2Yield  strength  at  -1100  K. 

^Density  at  -  room  temperature. 

-  Property  variation  over  temperature  range  of 
interest. 

^Tensile  strength  at  1100  K. 

^Nominal  value. 

^Above  300°C. 


containment  shell.  Based  on  the  Y-12  Plant  successful  experience  with 
304L  stainless  steel,  its  use  is  selected  for  testing  over  the  other 
stainless  steels.  Although  molybdenum  is  more  difficult  to  fabricate 
and  material  costs  are  much  higher,  it  appears  to  be  a  better  material 
than  stainless  steel  in  terms  of  strength  and  compatibility  with  lithium 
hydride.  The  compatibility  of  a  stainless  steel  shell  with  lithium 
hydride  nay  be  improved  with  a  thin  liner  of  pure  Iron  on  the  inside  and 
outside  shell  surfaces. 
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4.  EXPERIMENTAL  INVESTIGATION 

A  series  of  experiments  which  would  address  the  feasibility  of 
utilizing  encapsulated  lithium  hydride  as  an  energy  storage  medium  have 
been  identified.  The  series  of  tests  can  be  subdivided  into  five 
groups: 

1.  initial  scoping  experiments, 

2.  thin-shell  stress  tests, 

3.  thermal  characterization  tests, 

4.  hydrogen  diffusion  tests,  and 

5.  shell/lithium  hydride  compatibility  tests. 

To  date,  initial  scoping  experiments  have  been  performed  and  are  dis¬ 
cussed  in  Sect.  4.1.  The  other  four  groups  of  tests,  which  address 
specific  issues  discussed  previously,  have  not  been  conducted.  It 
appears  that  it  would  be  more  efficient  to  conduct  groups  of  experi¬ 
ments,  each  aimed  at  a  specific  issue,  rather  than  try  to  design  a 
"single"  experiment  to  answer  all  questions. 

Section  4.2  contains  a  discussion  of  experimental  and  analytical 
work  performed  to  investigate  the  feasibility  of  using  a  high  heat  rate 
graphite  induction  furnace  to  simulate  the  rapid  heatup  required  in  the 
actual  storage  system  ("real  system").  Rapid  heatup  of  capsules  in  the 
thin-shell  stress  tests  is  critical  and  is  highly  desirable  in  the  ther¬ 
mal  characterization  tests. 

4. 1  Scoping  Experiments 

With  the  assistance  of  Y-12  plant  personnel,  initial  scoping  exper¬ 
iments  were  performed  using  cylindrical  cans  containing  lithium  hydride. 
The  purposes  of  these  experiments  were  to  perform  preliminary  scoping 
tests  of  cylindrical  containments  and  to  gain  experience  in  the  thermal 
cycling  of  encapsulated  lithium  hydride.  The  tests  provided  insights 
into  container  (and  weld)  survivability  during  phase-change  and  the 
location  and  shape  of  the  void  formed  during  lithium  hydride  solidifica¬ 
tion.  In  addition,  experience  was  gained  in  the  areas  of  hydrogen  gas 


evolution  due  to  impurities  and  lithium  hydride  decomposition,  furnace 
equipment  functioning,  and  the  thermal  response  of  furnace  and  sample. 

The  experimental  setup  is  pictured  in  Fig.  4.1.  As  shown,  the  can 
fill-tube  is  attached  through  flexible  tubing  to  a  pressure  gage,  vacuum 
pump,  and  argon  purge  line.  The  lithium  hydride  container  was  placed  in 
a  stainless  steel  beaker  and  packed  with  MgO.  A  thermocouple  was 
located  adjacent  to  the  outside  of  the  lithium  hydride  container  about 
halfway  up  the  side.  The  beaker  was  placed  inside  the  resistance  heated 
furnace  (900  W  capacity)  and  the  top  of  the  furnace  was  covered  with  in¬ 
sulation.  The  cans  were  filled  with  preoutgassed  lithium  hydride  powder 
(chemical  analysis  of  similar  samples  indicated  99.16  mole  LiH). 
Following  filling  of  the  containers,  the  lid  with  fill-tube  was  welded 
to  the  can. 

Using  the  setup  described  above,  three  tests  were  performed.  Two 
tests  used  cylindrical  cans  2.54-cm  long,  3.81  cm  in  diameter,  and 
0.089-cm  (35-mil)  thick,  filled  with  about  14  grams  of  lithium  hydride. 
Based  on  density  considerations,  the  container  was  ~88%  full  at  the 
maximum  experimental  temperature  (~1050  K) .  In  the  first  test  (test 
C-l),  the  can  was  cycled  once  from  room  temperature  to  ~1000  K  and  then 
back  to  room  temperature.  In  the  second  test  (test  C-2),  the  can  was 
cycled  a  total  of  four  times  over  a  two  day  period  with  two  cycles  per¬ 
formed  each  day.  The  second  cycle  performed  each  day  was  initiated  from 
about  25  K  below  the  melt  point  rather  than  after  return  to  room  temper¬ 
ature.  However,  based  on  thermal  analysis  and  thermocouple  response,  it 
is  fairly  certain  that  complete  freezing  occurred  after  completing  the 
first  cycle.  Heating  from  about  25  K  below  the  melt  point  to  about  50  K 
above  the  melt  point  occurred  over  1  h;  the  cooling  process  was  also  1  h 
in  duration.  The  time-temperature  and  time-power  curves  for  both  tests 
are  shown  in  Figs.  4.2a  and  4.2b.  The  furnace  power  was  controlled 
(manually)  to  provide  heatup  from  ~25  K  below  the  melt  point  to  25  K 
above  the  melt  point;  a  similar  strategy  was  employed  on  cooldown. 

Since  the  powder  was  preoutgassed,  it  was  anticipated  that  only  a 
small  amount  of  gas  would  evolve  on  heatup  to  the  melt  temperature.  In 
all  tests,  this  was  observed.  In  heating  through  and  above  the  phase- 
change  temperature,  the  pressure  rose  in  tests  C-l  and  C-2  to  a  maximum 
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Fig.  4.3.  Calculated  equilibrium  hydrogen  pressure  curves  for 
est  C-2. 
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Fig.  4.4.  Void  and  lithium  hydride  cracks  formed  upon  freezing  in 
est  Cl. 


symmetrically  formed,  suggesting  that  the  can  was  not  level  during  test¬ 
ing  and/or  that  there  existed  some  nonuniform  heating/cooling.  Radio¬ 
graphs  of  the  test  C-2  can  show  a  similar  void  shape  and  location. 

Comparison  of  measurements  of  container  diameter  before  and  after 
C-l  and  C-2  showed  no  measurable  change.  However,  comparison  of  before 
and  after  top-to-top  measurements  for  test  C-2  showed  that  the  top-to- 
top  distance  had  decreased  about  0.05  cm  (20  mil)  (similar  data  was  not 
collected  for  test  C-l).  Since  the  top  and  bottom  of  the  can  were  flat, 
these  surfaces  could  easily  be  deformed,  possibly  by  "gripping"  or 
shrinkage  forces  generated  by  the  lithium  hydride  during  cooldown.  As  a 
result  of  heating  in  an  air  atmosphere,  cans  in  both  tests  oxidized. 
Based  on  visual  observation,  the  welds  performed  without  failure. 

Test  C-3  was  performed  with  a  cylindrical  can  with  the  same  dimen¬ 
sions  as  the  first  two  tests,  except  that  the  wall  thickness  was 
25  mils.  The  primary  objective  of  the  test  was  to  examine  void  behavior 
in  a  "hot-full"  container  (i.e.,  the  can  was  filled  with  16  grams  of 
lithium  hydride  to  produce  a  can  almost  completely  full  of  liquid  at  the 
highest  temperature  achieved).  Experience  with  lithium  fluoride  at 
TRWi&  and  ORNL57  indicates  that  in  some  cases  voids  tend  to  be  located 
near  the  center  of  the  container  following  solidification  even  in  the 
presence  of  Earth's  gravity.  Central  voids  occurred  when  the  containers 
were  both  "hot-full,"  and  cooled  from  the  entire  enclosing  surface.  TRW 
used  rectangularly  shaped  test  modules  (Fig.  4. 3)  to  study  lithium 
fluoride  phase-change  behavior.  Figure  A. 6a  shows  a  cross-sectional 
view  of  a  TRW  test  module  following  thermal  cycling  of  lithium  fluoride. 
Since  the  top  container  surface  had  a  significant  downward  dip  (caused 
by  machining  operations  in  making  thermocouple  well  grooves),  the 
lithium  fluoride  in  the  molten  state  contacted  the  top  surface.  Then, 
as  the  test  module  cooled,  lithium  fluoride  solidified  on  the  container 
surfaces  leaving  a  "centrally"  located  void.  Figure  4.6b  shows  a  con¬ 
trasting  situation  where  lithium  fluoride  did  not  contact  the  top  plate 
surface  in  the  molten.  In  this  case,  even  in  the  presence  of  fins,  the 
void  tended  to  form  at  the  top  of  the  container. 

Heat  storage  tubes  for  the  NASA  Brayton-cycle  heat  receiver  were 
filled  with  lithium  fluoride  at  ORNL.  A  photograph  of  a  tube  is  shown 


6a  source:  TRW  Power  Systems  Department,  "Brayton  Cycle  Cavity 
Receiver  Design  Study,"  TRW  Equipment  Laboratories, 
NASA  CR-54752  ER-6497,  Nov.  22,  1965,  Fig.  15,  p.  7 
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Fig.  4.7.  Bray ton-cyc le  heat  storage  tube. 

Fig.  4.7  source:  Gnadt,  P.  A.,  "Filling  Heat  Storage  Tubes  for  Solar 

Bray ton-Cyc le  Heat  Receiver  with  Lithium  Fluoride," 
Oak  Ridge  National  Laboratory,  ORNL/TM-27 32,  July 
1970,  Fig.  2,  p.  4. 
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test  C-3  was  slightly  lower  ('-80  mmHg).  Post-test  examination  revealed 
that  lithium  hydride  had  "crawled"  up  the  fill-tube  and  entered  the 
flexible  tubing  as  observed  in  test  C-2.  Although  not  apparent  during 
the  conduct  of  test  C-3,  lithium  hydride  leaked  from  the  can,  and. 
solidified  in  the  Mg"*.  The  MgO/kiH  mixture  formed  a  hard  material  which 
b  vided  tightly  to  the  test  can  and  beaker  making  it  verv  difficult  t.  n 
rv  :ove  the  can.  The  can  was  removed  hv  chiselling  the  ‘"gO  ’!  ill  mixture 
fr -m  the  car.  and  beaker. 

K;th  the  leakage  of  lithium  hydride  into  the  ‘V<>,  ami  the  1  •>*  .• 

l :  t  :;lun  hydride  tip  the  fill-tube,  .a  signi  f  i  cant  amount  -<f  U  t  h  * 

hydride  was  lost  from  the  can.  Radiographs  of  the  test  C-3  can  show  a 

v  >  i  .!  located  as  sketched  in  Fig.  4.9.  Apparently,  the  lithium  kvirilo 

■.as  able  to  "crawl"  up  the  can  walls  and  exit  through  the  open  i  op  ( s  1 

•.»: • .  re  the  weld  failure(s)  occurred. 

As  mentioned  previously,  in  tests  C- 1 ,  C-d,  tnd  C-  i  extensive 

..  ricking  was  observed  (visually  or  rad  iograph ica  1  ly )  in  the  cast  lithium 

' i y-griile  after  cooldown  (tor  example,  see  Fig.  4.4),  and  is  most  easilv 

; ::  r  ad  io  graphs .  hithiun  hydride  cracking  has  boon  observed  bv 

'th.-rs,  namely  Waldrop*  <  who  reports: 

1  ...  sore  cracks  are  always  present  in  f  in:  shed  f  U  tv> :  urn 
'■••irmie]  ;asf  ing  of  any  appreciable  sice  ...  we  have  found 
it  possible  to  prepare  cashing  of  considerable  size  which 
ire  sound,  in  the  sense  that  there  is  no  pipe,  although 
th.-re  .ire  numerous  cracks." 


'■ uiie;  by  Rapp  at  Thompson  Rano  Wooldrige,  Inc.,  as  reported  in 
In,  also  indicate  the  formation  of  lithium  hydride  cracks  follow- 


Fig.  4.9.  Sketch  of  void  shape,  size,  and  location  after  testing 
C-3  can. 


4.2 


Development  of  Large,  High  Heat  Pate  Furnace 


4.2.1  Graphite  induction  furnace  performance 

Comparison  between  the  capsule  experiments  performed  thus  far  and 
conditions  anticipated  in  the  "real  system"  show  some  verv  important 
differences,  other  than  the  obviously  snail  number  of  cycles  experimen¬ 
tally  performed.  As  indicated  previously,  operational  heatup  times 
from  '“700  to  1100  K  will  be  on  the  order  of  10  min,  whereas  the  lab¬ 
oratory  procedure  required  hours. 

In  order  to  obtain  "real  system"  heatup  tines,  an  existing  V—  1 2 
Plant  graphite  induction  furnace  will  be  modified  and  used  for  cycle 
testing.  The  furnace  is  cylindrical  with  a  peak  operational  power  of 
about  60  kW ,  and  about  0.01  m8  (0.4  ft8)  of  heated  volume.  The  size  and 
power  of  the  'urnace  will  permit  multiple  sample  testing. 

Initial  heatup  tests  were  performed  with  the  furnace  in  its  current 
configuration  to  assess  its  performance  capabilities.58  A  sketch  of  the 
furnace  cross-section  is  provided  in  Fig.  4.10.  As  shown  the  furnace 
contained  a  rather  large  graphite  sample  to  provide  some  thermal  mass. 
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Fig.  4.10.  Cross-sectional  sketch  of  cylindrical  graphite 
indu' “ i on  furnace. 


In  addition,  the  thermocouple  was  placed  inside  a  graphite  tube,  which 
also  provided  thermal  mass.  Figure  4.11  shows  the  thermocouple  response 
during  two  heatups,  and  the  following  cooldowns.  During  operation  the 
system  was  purged  with  argon  flowing  at  0.003  m^/min  (0.1  ftVmin). 
Following  a  brief,  low  power  pre-heat,  the  furnace  power  was  ramped 
quickly  to  60  kW.  In  the  first  cycle  (denoted  "Rl"  in  Fig.  4.11),  the 
thermocouple  temperature  rose  from  ~150°C  to  a  peak  temperature  of 
~915°C  over  about  10  min.  Furnace  power  was  shut-off  when  the  tem¬ 
perature  reached  850°C.  Following  cooldown  to  ~300°C,  a  second  cycle 
("R2")  was  performed  with  similar  results. 
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Fig.  4.11.  Graphite  induction  furnace  performance  in  current 
configuration. 

Fig.  4.11  source:  Wrenn,  George,  Y-12  Plant,  Martin  Marietta  Energy 

Systems,  Inc.,  Oak  Ridge,  TN,  January  1987. 
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These  data  indicate  tsat  rapid  heatups  can  be  obtained  in  the  in¬ 
duction  furnace.  It  is  anticipated  that  improvements  in  performance 
could  be  obtained  by  removing  some  of  the  thermal  insulation,  and  re¬ 
ducing  the  susceptor  wall  thicknesses. 

4.2.2  Projected  performance  of  furnace 

As  described  above,  checkout  tests  of  the  60  kW  induction  furnace 
indicate  that  a  graphite  sample  can  be  heated  in  the  furnace  at  a  rate 
of  about  100  K/min  (180°F/min).  In  order  to  estimate  the  thermal, 
response  of  a  sphere  containing  lithium  hydride  being  heated  in  the 
furnace,  a  computer  simulation  was  carried  out  assuming  the  furnace 
radiates  as  a  blackbody  initially  at  ambient  temperature  with  the 
furnace  temperature  increasing  at  a  rate  of  100  K/min  (180°F/min). 
Lithium  hydride  contained  in  a  7.62-cm  diameter,  0.0381-cm  (15-mil) 

thick  spherical  stainless  steel  shell  is  surrounded  by  a  muffle  which  is 
a  10.2-cm  inside  diameter,  0.076-cm  thick,  concentric  sphere  made  of 
niobium.  The  muffle  is  exposed  to  radiation  from  the  furnace.  All 
surfaces  are  assumed  to  have  an  emiss ivity  of  0.8,  and  all  heat  transfer 
is  by  radiation  with  no  heat  losses  from  the  muffle  or  lithium  hydride. 

The  resulting  temperature  profile  in  the  solid  lithium  hydride  at 
the  initiation  of  melting  is  shown  in  Fig.  4.12.  Also  shown  in 
Fig.  4,12  is  the  predicted  temperature  profile  at  the  initiation  of 
melting  resulting  from  suddenly  exposing  a  sphere  with  an  initial 
temperature  of  600  K  (620°F)  to  convection  heat  transfer  from  NaK  at 
1100  K  (1520°F),  and  the  corresponding  profile  obtained  in  a  low  flux 
furnace  (with  an  assumed  heatup  rate  of  200  K/h  or  360°F/h).  It  can  be 
seen  that  while  the  average  temperature  of  the  solid  in  the  high  flux 
furnace  sample  is  much  higher  than  that  expected  in  the  real  case,  much 
of  the  solid  is  cold  enough  to  retain  some  mechanical  strength.  Thus, 
some  of  the  stress-related  problems  previously  discussed  are  likely  to 
be  observed  in  tests  conducted  in  the  high-flux  furnace.  On  the  other 
hand,  in  low  flux  furnace  tests  the  solid  is  all  very  close  to  the  melt 
point  at  the  time  melting  begins  and  it  is  not  likely  that  stress 
problems  will  occur. 


5.  CONCLUSIONS 


Based  on  preliminary  system  concept  analysis,  encapsulated  lithium 
hydride  thermal  energy  storage  provides  signi f icant  benefits  to  space 
burst  power  systems.  The  use  of  encapsulated  I.i6H  in  a  packed  bed  with 
a  75%  packing  density  along  with  a  lithium  coolant,  will  provide  calcu¬ 
lated  system  mass  benefits  for  burst  periods  as  long  as  BOO  s.  In 
addition,  the  use  of  thermal  storage  will  permit  large  reductions  in 
radiator  area,  with  larger  benefits  at  shorter  burst  periods. 

A  group  of  feasibility  issues  associated  with  encapsulated  lithium 
hydride  thermal  energy  storage  have  been  identified  and  studied.  The 
feasibility  issues  include: 

1.  phase-change  induced  shell  stress  on  heatup, 

2.  hydrogen  diffusion  and  loss, 

3.  lithium  hydride  heat  transfer, 

4.  void  behavior/management,  and 

5.  material  considerations. 

The  key  issue  of  concern  is  the  possible  large  shell  stress  induced 
during  heatup.  Elastic  analysis  indicates  that  very  thick  shells  will 
be  required  to  prevent  shell  rupture,  assuming  the  lithium  hydride  forms 
a  structurally  sound  shell  following  solidification.  However,  cracks  in 
the  lithium  hydride  which  have  been  observed  to  form  during  cooldown, 
may  mitigate  shell  stresses  and  permit  the  use  of  a  thin  shell.  Never¬ 
theless,  analysis  has  shown  that  void  control,  via  sphere  heat  transfer 
control,  can  reduce  significantly  phase-change  induced  shell  stress.  In 
addition,  optimization  studies  have  indicated  the  benefit  of  high 
minimum  storage  temperatures  (^500  to  700  K),  at  which  lithium  hydride 
has  much  less  compressive  strength.  Alternatively,  a  flexible  shell, 
such  as  a  bellows,  may  be  viable. 

Based  on  material  consideration,  304L  stainless  steel  and  molyb¬ 
denum  are  leading  candidate  shell  materials.  Material  considerations 
included  material  compatibility  between  the  shell  and  lithium  hydride, 
lithium  and  hydrogen,  ductility,  strength,  density,  material  cost,  and 
ease  of  fabrication.  The  refactory  metals  niobium,  tantalum,  titanium, 
and  zirconium  have  been  eliminated  from  consideration  due  to  hydrogen 
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It  was  determined  that  SiC  is  not  compatible  with 
Lithium  at  —  1 L  00  K.  Other  ceramic  materials  have  been  eliminated 
primarily  due  to  their  poor  thermal  shock  resistance.  Molybdenum  and 
304L  stainless  steel  possess  reasonable  ductility  and  moderate  density. 
Molybdenum  has  greater  strength  and  better  material  compatibility  than 
304L  stainless  steel,  but  is  more  expensive  and  harder  to  fabricate. 

Based  on  projected  system  requirements,  a  0.0127-cm  (5-mil) 

stainless  steel  shell  will  provide  sufficient  hydrogen  containment  to 
prevent  significant  hydrogen  loss,  and  associated  loss  of  energy  storage 
density.  Due  to  molybdenum's  lower  hydrogen  permeability,  a  molybdenum 
shell  as  thin  as  1  mil  could  be  used  without  significant  hydrogen  loss. 

Initial  scoping  experiments  have  been  completed  for  cylindrical 
cans,  2. 54-cm  long,  3.81-cm  diameter,  and  0.0635-  and  0.0889-cm  (25-  and 
35-mils)  thick.  Four  thermal  cycles  were  successfully  completed  in  a 
low  heat  flux  furnace  for  a  0.0889-cm  (35-mil)  can.  Post-test  examina¬ 
tion  of  all  cans  tested  show  the  presence  of  numerous  cracks  in  the 
lithiun  hydride. 

Future  experimental  work  will  be  directed  at  determining  the 
feasibility  of  using  a  thin-wall  container.  To  this  end,  an  existing, 
high  power  graphite  induction  furnace  will  be  modified  and  used  to 
simulate,  as  closely  as  possible,  "real  system"  performance.  The 
material  and  fabrication  costs  associated  with  using  304L  stainless 
steel  and  molybdenum  will  be  investigated.  Presumably,  both  materials 
will  be  tested  and  will  permit  not  only  an  assessment  of  shell  heatup 
stress,  but  also  hydrogen  loss,  material  compatibility,  and  lithium 
hydride  heat  transfer. 

Future  analytical  work  will  extend  current  capability  to  include 
the  effect  of  surface  tension  gradients  on  void  shape,  void  location, 
and  liquid  flow,  and  boundary  layer  effects  on  the  void.  In  addition, 
the  void  and  natural  convection  models  will  be  integrated  into  the  beat 
transfer  code.  These  models,  in  addition  to  nodels  for  stress  analysis 
and  hydrogen  loss,  can  then  be  compared  to  data  to  determine  consis¬ 
tency. 

Since  the  micro-gravity  space  environment  will  produce  different 
behaviors  than  those  on  Earth,  some  basic  verification  tests  should  be 


performed  in  space  to  determine  the  accuracy  of  model  assumptions  and 
results,  and  whether  all  important  differences  have  been  anticipated. 
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